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Abstract—Equations are developed for calculating the thermodynamic efficiency of a chemical 
manufacturing process, using as a basis for the efficiency the second law, rather than a simple 
energy balance. This requires a consideration of the temperature at which heat is recovered, as 
well as the amount recovered and involves the use of the “ availability * function, as developed 
by Keenan. The thermodynamic efficiency of a typical ammonia oxidation process is shown 
to be about 6%, or is 11% when prescribed conditions of irreversibility are allowed for. By 


wm contrast the ordinary thermal efficiency, based on the first law, is as high as 32%. 
/57 Résumé—L’auteur développe des équations pour le calcul de lefficacité thermodynamique 
956 d'un procédé de fabrication chimique en utilisant le second principe plut6t qu'un bilan thermique. 


Ceci nécessite la connaissance de la température 4 laquelle la chaleur est récupérée aussi bien 
que la quantité récupérée et implique l'utilisation de la fonction “ disponibilité ” développée 
par KEENAN. 

L’efficacité thermodynamique d'une oxydation classique de l'ammoniac est d’environ 6% 
et de 11% quand il est tenu compte de l'irréversibilité. Au contraire, l'efficacité thermique 


basée sur le premier principe est de ordre de 32%. 


Waste Worx 


Tue notion of the efficiency of a heat engine is 
a very familiar one ; the second law sets a certain 
limit on the conversion of heat into work and 
the efficiency of any irreversible engine may be 
calculated relative to the reversible engine 
considered as a “ target.”” The corresponding 
notion of the thermodynamic efficiency of a 
process for producing chemicals is much less 
familiar and this is the subject of the present 
paper. We shall be concerned with the target 
against which any actual process should be 
assessed, using as our basis the second law 
rather than a mere energy balance. 

Let it be supposed that the temperature, 
pressure and composition of the raw materials 
have already been decided upon, and also the 
temperature, pressure and composition of the 
required product. The initial and final states 
of the process materials having been fixed in 


this way, their changes of internal energy, 
entropy and volume are completely determined. 
These changes will be denoted AU, AS, and 
AV respectively. 

As a measure of the efficiency of this process 
we can use either the minimum input of work, 
or the minimum input of heat from a source at 
a specified temperature. The former is the 
more convenient and we shall therefore use, 
as our criterion the least amount of mechanical 
energy which must be supplied when, at the 
conclusion of the process, the only systems 
which have undergone any change are the process 
materials and the mediumt. By the latter is 
meant the atmosphere and the earth’s surface, 
(or cooling water at the temperature of the 
earth). 

That is to say, the efficiency is to be computed 
on the basis that, at the end of the process, 
all auxiliary heat reservoirs, such as quantities 


*Based on two lectures given in the Chemical Engineering Department of the University of Minnesota in 1953. 
tAnd, of course, the position of the “ weight" which measures the performance of the mechanical energy. 
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of steam, have been restored to their original 
states, ie. the medium is to be taken as the 
only source or sink. This is clearly the only 
sound basis, since the medium is the only heat 
reservoir available to man which does not require, 
for its maintenance, the using up of some of his 
resources, 

In a considerable number of chemical processes 
it is possible, at least in principle, to obtain an 
output of useful energy. In such cases the 
corresponding criterion is the marimum amount 
of mechanical or electrical energy which could 
have been obtained from the specified process, 
when the only other body which undergoes 
any change is the medium. Of course, in the 
special case where the initial and final states 
of the process materials are at the same tempera- 
ture and pressure as the atmosphere, the minimum 
work which must be applied is equal to the 
increase of their Gibbs free energy. (Or the 
maximum work which might be obtained is 
equal to the decrease), But in cases where 
the temperatures and pressures are not the same 
as those of the atmosphere the minimum work 
will be shown to be equal to the increase in 
* availability " of these materials. 

The thermodynamic function known as the 
availability was developed by Kren an [1] on the 
basis of some work of WiLLarp but has 
not been applied at all extensively to chemical 
processes. The significance of this function for 
the present purposes will now be outlined. 

Consider first of all the process as it is actually 
carried out. The first law may be written 


AU =q — x’, 


where AU is the increase in internal energy 
of a fixed quantity of process materials passing 
through the chemical plant, w’ is the work done 
by these materials and q is the heat they take in. 
For the present purposes it is convenient to 
write this equation in more detail and to change 
the sign of the work because most real processes 
require a positive input [e.g. to pumps]. 
Thus 

AU = + Po AV +e, (1) 
where 
Yo is the heat given to the medium ; 
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2q, is the sum of all other heat effects 
(e.g. from various supplies of steam), 
taken as positive for heat taken 
in by the process materials ; 

P,AV is the work done in displacing 
the atmosphere at its constant 
pressure Py; 

w is the total of all other items of 
work, taken as positive for work 
done on the process materials. 

This equation will apply whether the process 
is batchwise or continuous, provided that w 
is taken to include the work required to cause 
inflow and outflow and is not merely the shaft 
work within the process itself. (Otherwise AU 
would be replaced by AH in the case of steady 
flow. 

According to the second law the sum of all 
entropy changes caused by the process can only 
be positive or zero. Thus 
AS + AS, + 2 AS; > 0, 

AS is the entropy change of the process 
materials ; 

AS, is the entropy change of the medium; 
SAS, is the sum of the entropy 
changes of the reservoirs which provide 
the heat effects q;. 


where 


Alternatively it is convenient to define a 
quantity o by means of the identity 


AS + AS, + Z AS, =o, (2) 


so that o means literally the created entropy. 
Then according to the second law 


o> 0. (3) 


This device of writing a symbol for the amount 
by which the entropy of the universe has increased, 
as a result of the process, was first used at all 
extensively by the Belgian physicist pe DonpeEr. 
It introduces nothing which is not already 
contained in the second law but it does assist 
in making problems of irreversibility more 
easy to survey. 

Now the medium acts merely as a heat reservoir 
at constant temperature 7, (which is in the 
region 280-300°K). Hence 


AS, = qo/To 
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and therefore (2) may be written 
AS + FAS, (4) 


Eliminating g, between equation (1) and (4) we 
obtain 


AS + (- AU + 24, +w — P, AV) 
+ ZAS,;=0 


Multiplying through by 7, and rearranging 


w+ (9; + T, AS,) = + AU + P, AV 
— 48 (5) 


So much for the process as it is actually carried 
out. Consider now what would be the work 
requirement if no heat reservoir, other than the 
medium, were allowed to undergo any change. 
The terms g; and AS, in equation (5) would 
therefore be zero. On the other hand, the last 
three terms on the right hand side of (5) will 
have precisely the same values as for the actual 
process, since the change in state of the process 
materials is taken as fixed. For the same amount 
of created entropy the total work input w, of 
such a process would therefore be given by the 
right hand side of (5). Thus 


w,= + AU + AV — T, AS (6) 


The quantity w + 2(q; + T, AS;) on the left 
hand side of (5), which is the same as w, in (6), 
is therefore the total amount of work which 
would be needed for the given process, for a 
particular value of o, the created entropy, when 
the only systems which undergo any change 
are the process materials and the medium. 
This may perhaps be seen more clearly by noting 
that, in the process as actually carried out, 
2 (q; + Ty AS,) is the additional work needed 
to restore the various heat reservoirs to their 
original states, at the expense of the medium. 
For consider the restoration of heat dq; to the 
ith reservoir, by use of a refrigeration cycle, using 
the medium as the source. The necessary work 
is given by the familiar Carnot expression : 
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where 7; is the momentary temperature of the 
reservoir. Thus the work needed to restore 
q, to the reservoir is 


=q + T, AS; 


(where the sign of AS; arises from the fact that 
AS, refers to the process of taking heat from 
the reservoir). 

It follows that the quantity w + I(q; + T, 
AS;), which is denoted w, in (6), may be referred 
to as the “ total equivalent work requirement” of 
the process, for a given value of ¢. Now according 
to the second law o can only be positive or zero. 
Therefore the least work requirement of the 
process is given by 


= AU +P, AV—T,AS, (7) 


and this occurs when the process is carried out 
reversibly. 

Conversely if the signs of AU, AV, and AS 
are such the process is able to deliver mechanical 
work, then the above expression gives the 
negative of the maximum work which may be 
obtained, the only changes being in the process 
materials and in the medium. 

Comparing equations (5) and (7) [or (6) and (7)] 
it is evident that the quantity T,o is the 
additional quantity of work* which is used in the 
actual process as compared to the ideal reversible 
process. Thus Tyo may be called the * wasted 
work ” or * dissipated energy.” It is identical 
with the quantity which Krenan calls the 
“irreversibility” of the process. For the present 
purpose the term waste work is the more descrip- 
tive. One of the functions of the chemical 
engineer is to device his process so that o is as 
small as possible, consistent with not too large a 
capital expenditure on the plant. 

The quantity on the right hand side of equation 
(7) has been called by Keenan the increase in 
availability in the process in question. DopcrE 
[3] has given it the symbol AB, 


AB = AU + P, AV — T, AS. (8) 


*Or the additional quantity of heat, expressed as its 
work equivalent by use of the Carnot factor. 
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Of course, if the initial and final temperatures 
and pressures of the process materials were 
the same as those of the medium, AB would 
be the same as AG, the increase of Gibbs free 
energy and we should have the familiar result 
that the minimum work which is needed is 
equal to the increase of G. However the utility 
of the present equations is not limited to isother- 
mal or isopiestic processes, (as is the case when 
the Gibbs free energy itself is used). The equations 
may be employed, for example, to calculate the 
minimum necessary input of work in cases where 
the reagents enter the plant at a temperature or 
pressure quite different from those of the out- 
going products. 

Whenever the product is a final product, 
in the sense that it is not to be used in some 
subsequent process, the most economic conditions, 
from a thermodynamic standpoint, are obtained 
by discharging this product at the same tempera- 
ture and pressure as the medium. For if it 
were at a temperature T different from Ty, an 
output of work could be obtained by operating a 
heat engine between the product and the medium, 
using whichever is the hotter as the * source ” 
and whichever is the cooler as the “ sink.” 
Similarly if the product were discharged at a 
pressure different from Py, the pressure of the 
atmosphere, additional work could be obtained 
by operating an expansion engine. 

To summarise, equations (5) or (6) may be 
written 

w,= Tyo + SB (9) 
or in words 


(total equivalent work requirement of process) 
== (wasted work) + (increase in availability 
of process materials). 


This increase in availability, AB, is entirely 
determined by the initial and final states of the 
materials, together with the temperature and 
pressure of the medium. However the engineer 
still retains the liberty to try to make o, and 
therefore w, as small as possible. This requires 
a suitable choice of the path between the given 
initial and final states, ie. the choice of the 
details of operation. Moreover it is evident 
that in processes for which AB is negative it is 


possible to obtain an output of work, That is to 
say w, is negative if AB is negative and if o 
is small or zero*. 

It is evident that the smallness of To, 
relative to AB, is a measure of the efficiency of 
the process. However it is not intended to 
propose an actual definition of the efficiency, 
since it seems impossible to put forward a defini- 
tion which is equally satisfactory for AB positive 
or negative. In the example which is now to 
be discussed it will be clear enough how the 
efliciency is to be understood. 


2, Tue Waste Work or AMMONIA 
OXIDATION 


As an illustration of the preceding discussion 
it is instructive to consider the production of 
nitric acid from ammonia. This is a process for 
which AB is negative and the process is capable 
of yielding a far larger output of useful work 
than is actually obtained. 

It will be supposed that the ammonia is 
delivered to the process as a liquid at 20°C, 
under its vapour pressure of 8-5 atm. The other 
reagents are water and air, each at 15°C and 
latm. The end product of the process is 60% 
(wt.) nitric acid, together with the unused 
components of the air, all of which are assumed 
to be at 25°Ct and one atmosphere. 

For simplicity it will be supposed that the 
chemical efliciency of the conversion of ammonia 
into nitric acid is 100%. That is to say, for the 
purpose of an analysis of the thermodynamic 
efliciency of the process we shall disregard the 
complicating effect of the chemical efficiency, 
which is usually rather less than 100°%,**. 


*As a corollary, a spontaneous process within the 
medium is one in which there is a decrease of availability. 


+The various temperatures are chosen as slightly 
different from each other in order to illustrate the method 
of calculation for cases where the differences are much 
larger. 

**The main chemical inefliciency is the conversion of 
about 5% of the ammonia into nitrogen, instead of into 
nitric oxide, on the platinum catalyst. Allowing for this 
would not appreciably alter the second law analysis of 
the process as given in this paper. 
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The overall chemical process may therefore 
be written 


1-34 H,O + NH, (1) + 9-6 air 
= HNO, + 2:34H,O +76N,, (A) 


where the 1-34 moles of water on the left hand 
side of the equation is that which is required 
to give a 60%, acid as the final product on the 
right hand side. The 7-6 moles of nitrogen 
on the right-hand side is, of course, that quantity 
of nitrogen which is incidental to the usage 
of 2 moles of oxygen in the reaction process. 
Standard values at 25°C of the free energies 
and enthalpies of formation of the various 
compounds which occur in the above equation 
are available in the literaturet. Using these 
data, together with the known partial pressures 
of H,O and HNO, above the 60%, acid (in order 
to allow for their activity coefficients), it may be 
calculated that the following are the changes 
of the Gibbs free energy and enthalpy, when one 
gram mole of nitric acid, as a 60%, solution, is 
formed at 25°C from air and from liquid ammonia. 
under its vapour pressure, at this temperature, 


AGog, = — 742 kg cal 
AH, = — 98-2 kg cal 


The volume change in this process is almost 
exactly that which corresponds to the disap- 
pearance of 2 moles of oxygen at 25°C and one 
atmosphere. 
Hence 
A (PV = — 48°8 1 atm. 
= — 1-2kgcal 
Therefore 
AU gg, = AH — A(PV) 
= — 97-0 kg cal 
ASoog = (AH — AG)/T 
= — 0-0805 kg cal per °C. 


These are the changes in AU and AS which 
pertain to an isothermal process at 25°C. The 
actual process to be discussed was stated above 
to be one in which the products are at this 
temperature, whilst the liquid ammonia is 


+The values used are from the National Bureau of 
Standards compilation, 1952. 


under its vapour pressure at 20°C, and the air, 
together with the 1-34 moles of water needed 
to give the 60%, acid, are at 15°C. To the above 
values of AU and AS we must therefore add 
the changes of U and S which are associated 
with the heating of liquid NH,, under its vapour 


pressure, from 20°C to 25°C and the heating of 
9-6 moles of air and 1:34 moles of water from 
15°C to 25°C. 

When these small corrections are made, we 
finally obtain the following values for the changes 
of internal energy and entropy, between the 
assigned initial and final states of the process 


AU = — 96-3 kg cal 
AS = — 0-0792 kg cal per °C. 
Also for the volume change 
AV = — 41 1. 


Therefore the change of the availability is 


given by 
AB = AU + P, AV — T, AS 


— 96-3 — 1-0 + 22-8 
= — 74:5 kg cal, 


where P, and T, refer to the medium, and have 
been put equal to one atmosphere and 288°K 
respectively*. 

It follows that the process is capable of pro- 
viding 745kgcal of mechanical energy (for 
example as electrical power) per gram mole of 
nitric acid produced. This requires reversible 
conditions of operation, and a good approximation 
to this might be attained by carrying out the 
reaction electro-chemically, in a galvanic cell, 
if this were practicable. 

Consider now the process as it is actually 
carried out, between the assigned initial and 
final states. Some typical figures, which refer 
to the customary process in which ammonia is 
oxidised on a platinum-rhodium catalyst, followed 
by absorption of the nitrogen oxides in packed 


*In this example the difference between AG and AB 
is seen to be very small. Much larger differences are to be 
expected in examples where the process materials enter 
or leave the system at temperatures appreciably different 
from Tp. 


towers at atmospheric pressure, are as follows* 


69 kwh per short 
ton HNO,. 


0-8 tons per ton 
HNOs. 


The second figure refers to the output of steam 
from a waste heat boiler which takes the hot 
gas from the ammonia oxidation. The steam in 
question is saturated at 50 lb. per square inch 
and is obtained from feed water at 20°C. 

These figures may be re-expressed as follows 


input of electrical power 


steam recovery 


electrical power usage 4-1 kg. cal per 


g mole HNO, 
31-9 kg cal per 
g mole HNO, 


heat transferred to steam 


Now it is to be noted that the energy which is 
recovered from the process, 31-9 kg cal, although 
quite substantial, is recovered as heat and not 
as work. If the local conditions at the factory 
are such that there is a requirement for heat 
(e.g. for space heating) this amount of energy 
could, in a sense, be allowed to count at its full 
value. However, when we are examining the 
absolute efliciency of the process, in the meaning 
of the second law, the fact remains that we 
could have obtained work from the process and 
actually we have only obtained heat — and work 
is always more valuable than heat on account of 
the Carnot conversion factor. 

In fact the work equivalent of the recovered 
heat is the negative of the term (g; + Ty AS;) 
on the left hand side of equation (5)t. Here q; 
has the value — 31-9kgceal (negative since 
heat has been given to the steam in the process 
considered) and AS; is the increase of entropy 
in the process : 


(feed water at 20°C) + (saturated steam at 
50 Ib. per sq. inch). 


Since there are 0-8 « 63 = 50-4 grams of steam 


*These figures were kindly supplied by the Directorate 
of Ordnance Factories of the Ministry of Supply, but are 
entirely typical of ammonia oxidation as it is normally 
carried out. 

+Note that there is only one such term in the process 
considered. 
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per gram mole of nitric acid, this entropy change 
is readily calculated from steam tables and, 
multiplying by Ty. we obtain Ty, AS; = 
23-1 kg cal. Therefore, 


work equivalent of recovered heat = 
— (— 31-9 + 23-1) 
= 88 kg cal per mole HNO,. 


We now have all the terms occuring in equation 
(5) except Tye and this is obtained by difference : 


w + T, AS,) = T,0 + AB 
+1 — 8-8 = T,o — 745 


Therefore T, o has the value 69-8 kg cal per mole 
of HNO,. 

In brief there is a very large waste work or 
dissipation of energy in this particular process. 
The work equivalent of the recovered heat, less 
the input of electrical energy, is 88 — 41 = 
47 kg cal; a perfectly reversible process would 
have yielded as much as 745 kg cal of useful 
work. The efficiency of the process, in the sense 
of the second law*, might therefore be defined 
as the ratio of these two figures and is 6%. 


3. Prescripeo Decree or IRREVERSIBILITY 


The efficiency of the process, as computed 
above, could be raised to 100°, only if some 
means were available for carrying out the reaction 
under reversible conditions (for example in a 
galvanic cell using ammonia and oxygen elec- 
trodes). For this reason the efficiency of 6% 
may be called the absolute efliciency, since it is 
calculated relative to a standard of perfect 
reversibility. 

On the other hand it is pertinent to ask the 
question : if, for practicai reasons, it is prescribed 
that the reaction must be carried out in a certain 
way-—namely by irreversible oxidation on 
platinum at 800-900°C — what is the efficiency 
with which the actual process is carried out, 
within the framework of this limitation ? 


*Not to be confused with the efliciency in the sense of 
the first law, as based on the energy balance. This would 
normally be defined as the ratio of the heat which is 
actually recovered to that which could have been recovered 
if there were no loss by convection and radiation to the 
atmosphere. This efliciency is 32°, in the above example. 
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A similar situation exists in a great many 
industrial reactions and it is useful to calculate a 
practical efficiency, in addition to the absolute 
efficiency discussed above. The practical elli- 
ciency may be defined in such a way that it is 
unity (or 100%) when certain parts of the 
process are carried out under prescribed con- 
ditions of irreversibility, and all other parts are 
carried out under ideal reversible conditions. 

As an example we shall consider once again 
the oxidation of ammonia to give 60°, acid, the 
initial and final states being as described pre- 
viously. The practical efliciency will be defined 
in such a way that it is 100°, when the following 
conditions are satisfied. 


(a) All parts of the process are carried out 
reversibly except the chemical reactions 
and the mixing of the reagents. 


(b) The first reaction, the catalytic oxidation 
of ammonia, takes place at 850°C 
(1,123°K) and 1 atm. pressure and is 
carried out under completely irreversible 
conditions in the sense of performing 
no work (other than the displacement of 
the atmosphere). Similarly the second 
reaction, the absorption of the oxides of 
nitrogen, will be specified as taking place 
at 25°C and 1 atm. Under these con- 
ditions the heat evolved in the reactions 
is given by the change of a function of 
state and is — AH*. 


(c) The heat of reaction is recovered as 
mechanical work at the Carnot conver- 
sion factor (i.e. in a reversible engine). 


*It is important that the heat evolution in a reaction is 
determined by the change of a function of state of the 
process materials only when one or the other of the 
following conditions is satisfied. (a) The process is 
reversible and isothermal. In this case the heat evolution 
is — T AS (and is not — AH). (b) The process is irrever- 
sible and the only work done is a displacement of the 
environment at constant pressure (a condition which 
prescribes the amount of work). In this case the heat 
evolved is again determined by the initial and final 
states and is — AH. (c) The process is irreversible and 
takes place at constant volume (zero work). The heat 
evolved is — AU. 
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The amount by which the practical elticiency 
is less than 100%, is therefore a measure of the 
departure from reversibility of parts of the 
process other than the chemical reactions, e.g. 
the pumping of gases and acids, conversion of 
heat of reaction into useful work, etc. 

The practical efficiency may be calculated by 
either of two methods— 


(1) A calculation is carried out on the amount 
of work which may be recovered from 
each part of the process, subject to the 
above conditions concerning irreversibility 
in the reactions. 


(2) Instead of calculating the recovery of 
work, we calculate how much is “ wasted,” 
i.e. the value of Tyo in the irreversible 
parts of the process. 


Although the results are equivalent, an under- 
standing of a given manufacturing process can 
best be attained if both methods are used. 
The first shows where to look for sources of 
recoverable energy, whilst the second emphasises 
the causes of wastage. Since AB for the overall 
process is known, the two answers can be cross- 
checked and any discrepancy will indicate that 
some factor has been overlooked. 

The first type of calculation is very familiar 
and therefore only the second will be discussed 
here. For this purpose it is first of all necessary 
to show how the entropy creation in a reaction 
may be calculated. 

Consider any reaction taking place at constant 
temperature and pressure, the heat of reaction 
being transmitted to some other body (e.g. steam) 
which is at the same temperature. The created 
entropy is 

o = AS + AS,, 


where AS and AS, refer to the entropy increases 
of the reaction system and the other body, 
respectively. If the only work done by the 
reaction system is due to its change of volume at 
constant pressure, then the heat transmitted to 
the other body is equal to — AH, the decrease 
of enthalpy of the reaction system. The last 
equation may therefore be written 
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(10) 


The entropy created in a chemical reaction* 
taking place at constant temperature and pressure 
is therefore equal to — AG, the decrease of Gibbs 
free energy of the reacting system, divided by 
the temperature T at which feaction takes place. 
The “wasted work,” or “ dissipated energy,” 
or “ irreversibility,” of the process is therefore 


_ 1, 4G 


= T 


(11) 


4. Exampte on Ammonia OXIDATION 
The overall reaction (A), as written in Section 2, 
must now be separated into its two component 
parts, which take place on the platinum catalyst, 
at 850°C and 1 atm, and in the absorption 
towers, at 25°C and 1 atm., respectively. These 
aret 


NH, + 96.air = NO +5 H,0(g) +0, 


+76N, (B) 


1-84 H,0 (I) +NO + 5H,0(g) 
+ 50s +7-6N, = HNO, 


+ 2-84H,0 + 7-6N,, (C) 


where the 3/2 moles of gaseous water on the right 
hand side of (B) are formed in the oxidation 
reaction and the additional 1-34 moles of liquid 
water on the left hand side of (C) is the additional 
water which must be added to form the 60°, acid. 

From data in the literature on the free energy 
and enthalpy change in reaction (B) at 25°C, 
together with heat capacity data, it may be 
calculated that the value of AG at 850°C (1,123°K) 


*For methods of calculating the entropy production 

in other types of natural process see the author's 
Thermodynamics of the Steady State, Methuen Monograph, 
1951. 
+ The writing of the equation in this manner, which in- 
cludes inert components such as nitrogen, assists in the 
calculation of the correct free energy charge, which in- 
cludes a small term due to the free energy of mixing. 


is about — 78 kg cal/g mole NO. Using equation 
(11) the waste work is therefore 


288 
= —— 78 
7123 


= 20 kg cal /g mole NO. 
Similarly in reaction (C) 


288 


998 


as in equation below 
= 12-3 kg cal/g mole HNO . 


Adding together the two items of waste work 
gives a total of 32-3 kg cal per mole of HNO, and 
this is therefore the amount of energy which 
inevitably we fail to recover as useful work when 
the two reactions are carried out under the 
prescribed conditions of irreversibility. 

If the whole process had been carried out 
reversibly it would have been possible to obtain 
74-5 kg cals. of useful work, this being the value 
of — AB for the overall process, as calculated 
previously. Therefore the maximum amount of 
mechanical energy which could possibly be 
obtained under the prescribed conditions is 


74:5 — 32-3 = 42-2 kg cal/g mole HNO, 


(as could also have been calculated by Method 1 
on page 7). 

According to the figures quoted in Section 2, 
the industrial process, as it is actually carried 
out,’ uses an input of 4:1 kg cals. of electrical 
energy (for pumping), whilst the recovered 
steam has a work equivalent of 8-8 kgcal. The 
net recovery is therefore 


8-8 — 41 = 47 kg cal. 


The practical effiiciency of the process may 
therefore be expressed by the fraction 4-7/42-2 
or 11%. 

This figure is considerably higher than the 
absolute efficiency of 6%, as calculated in 
Section 2, but is still very low. This is due to 
several factors. 
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(a) The need for an input of 4-1 kg cals. of 
electrical energy, for the purpose of 
forcing the gases through the system and 
also for acid pumping. In brief there is an 
irreversibility due to fluid friction which 
could be avoided only by eliminating the 
pressure drop. 

(b) The incomplete recovery of the heat of 
reaction, the remainder being lost to the 
atmosphere and to the cooling water. 


(c) The fact that the steam is obtained at a 
temperature of only 138°C* and not at 
the highest possible temperature of 850°C. 
The latter is feasible, in principle, in view 
of the fact that the ammonia oxidation 
takes place at this temperature. 

This latter factor greatly reduces the work 
value of the steam, on account of the Carnot 
conversion factor as discussed previously. 

Of course if the efficiency were computed on 
the basis that the steam is to be used for heating 


*Saturated stream at 50lb. per square inch. 


purposes, the situation would appear considerably 
better. However this would be equivalent to 
computing the efficiency on the basis of the first 
law, rather than of the second. 

In conclusion it is perhaps useful to summarize 
the main essentials for keeping entropy creation 
at a low level. 


(a) All heat transfers should take place at the 

least possible temperature difference. 
(b) Pressure drop and other forms of friction 
should be kept at a minimum. 


(c) Chemical reactions should be carried 
out under “ resisted ’’ conditions, so that 
they yield useful work. 


The latter factor is certainly one of the most 
important in the chemical industry but unfor- 
tunately it is a matter of the greatest practical 
difficulty to carry out reactions under approxi- 
mately reversible conditions. If some of the 


large-scale industrial reactions could be set up 
in the form of a fuel cell it would undoubtedly 
lead to a considerable saving of energy. 
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Abstract— It is well known that when petroleum fractions of different boiling range are mixed, 
their volumes are not strictly additive. However, little effort has been made to investigate 
the effect of composition, chemical nature of the components, and temperature on the magnitude 
and sign of the deviation from additivity. In order to obtain information on these factors, but to 
avoid the difficulties attendant with the interpretation of the data on such a complex mixture 
as petroleum, a study was made of the ternary system, cyclohexane-n-heptane-benzene. Volume 
changes on mixing were computed from the specific volumes of the pure compounds and their 
binary and ternary mixtures which were determined experimentally at 60°F and 100°F. 

An increase in volume on mixing was noted in all mixtures studied, with the value of the 
maxima at 60°F ranging from about 0-30%, in the case of a 60-40 mixture (by weight) of cyclo- 
hexane and n-heptane to 0-70 %, in the case of a 50-50 mixture of cyclohexane and benzene. In the 
ternary system a maximum of about 0-7% increase in volume was found in mixtures containing 
equal proportions of cyclohexane and benzene and up to 20 weight %, n-heptane. The magnitude 
of the effect was reduced by an increase in temperature. It was found that the experimental 
results could be represented satisfactorily by equations based on the theory of regular solutions. 

The data are of interest in connection with the manufacture of certain petroleum products 
that are prepared by the blending of different fractions. 


Résumé— Dans un mélange de fractions de pétrole de points d’ébullition différents, le volume 
résultant n‘obéit pas a la loi d'additivité. Quelques recherches, cependant, ont été faites pour 
étudier Veffet : de la composition, de la nature chimique des composants, de la température 
sur importance et le signe de lécart a ladditivité. Pour obtenir des renseignements sur 
ces facteurs et dans le but d'éviter les difficultés inhérentes 4 linterprétation des données sur 
un mélange aussi complexe que le pétrole, les auteurs ont étudié le systéme ternaire eyclohexane- 
n-heptane-benzéne. Les variations de volume du mélange ont été estimées a partir des volumes 
spécifiques des constituants et des meélanges binaires et ternaires expérimentaux a 60 et 100 °F. 

Dans tous les mélaunges étudiés, les auteurs ont remarqué une augmentation de volume 
maxima & 60 F et variant de 0-3°,, pour un mélange de 60-40 (en poids) de ceyclohexane-n-heptane 
& 0-70% pour un mélange de 50-50 cyclohexane-benzéne. Dans le systéme ternaire, augmenta- 
tion de volume est de 7°, dans des mélanges cyclohexane-benzéne en égales proportions 
contenant jusqu’éA 20% en poids de n-heptane. Un accroissement de température réduit 
Vimportance de l'augmentation. Les résultats expérimentaux peuvent se représenter d'une 
maniére satisfaisante par des équations basées sur la théorie des solutions réelles. 

Toutes ces données présentent un intérét pour la fabrication de certains produits pétroliers 
préparés par mélange de différentes fractions. 


Ir 1s common knowledge that when two dissimilar is insignificant. While this may be true in the 
liquids are mixed a change in volume and in case of certain specific blends it is not generally 
temperature occurs. In general, the more dissimi- true. TureLe and Kay [13] have shown that in 
lar the compounds the larger will be the thermal the blending of a 50-50 mixture by weight, of 
and volumetric effects. In the blending of petrol- light and heavy naphtha, a volume shrinkage of 
eum stocks, the similarity in the nature of the about 0-22% occurs. Iorre [5] determined the 
component hydrocarbons has led to the tacit change in volume on mixing in nine binary hydro- 
assumption that the volume change due to mixing carbon systems and in five petroleum systems 
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each formed by blending various amounts of two 
different petroleum fractions. He reported volume 
changes amounting to as much as 1%, when the 
fractions blended, differed widely in their mole- 
cular weights and physical properties. ConNoLLy 
[2], in a study of the n-butane-benzene-cyclo- 
hexane systems, reported volume shrinkages in 
50-50 mixtures of the n-butane-benzene and 
n-butane-cyclohexane binaries that varied from 
0-4% at 60°F to 15% at 140°F. Reeves [7] and 
Kurtz and Sank1n [6], in summaries of the data 
on the volume change on mixing of hydrocarbon 
liquids reported in the literature, have presented 
additional evidence that volume changes on 
mixing in hydrocarbon systems can be appreciable 
and, in some blending operations, may affect the 
economic values assigned to individual products. 

To obtain further information on the sign and 
the magnitude of the volume change on mixing 
in hydrocarbon systems, and on the factors 
affecting the change, a study of the ternary system 
cyclohexane-n-heptane-benzene was made. This 
system which is composed of compounds repre- 
sentative of the three general classes of hydro- 
carbons occurring in petroleum fractions, was 
chosen with the objective in mind of determining 
principally, the effect due to differences in the 
chemical nature of the components. Volume 
changes on mixing were computed from the 
specific volumes of the pure compounds and their 
binary and ternary mixtures which were deter- 
mined in glass pycnometers at 60°F and 100°F. 
Equations for representing the relation between 
the percentage volume changes and the composi- 
tion at constant temperature in binary and ternary 
mixtures were obtained by the application of the 
theory of regular solutions developed by HiLpE- 
BRAND [8] and Scatcnarp [9], [10], and the 
calculated values were compared with the experi- 
mentally derived values. 

Volume changes in the binary system, cyclo- 
hexane-benzene have previously been determined 
by Scatcuarp, Woop, and Mocue [11] at 80°C 
and later by Woop and Austin [14] from 15°C 
to 75°C at 10°C intervals. Brown and Ewatp 
[1] have obtained similar data for the n-heptane- 
benzene system at 25°C. Kurtz and Sanxin [6] 
have reported volume change for one mixture of 
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cyclohexane and n-heptane. So far as is known, 
ternary mixtures of cyclohexane, n-heptane, and 
benzene have not been studied before. 


EXPERIMENTAL 


For the determination of the specific volume, three 
pycnometers were used so that simultaneous determina- 
tions on three samples could be carried out. Each of the 
pycnometers was constructed from a 50 ml pyrex volum- 
etric flask provided with a ground glass stopper, by 
replacing a section of the neck with a short piece of 
2mm bore capillary. A ring was etched around the 
capillary at about its mid-section to serve as the mark 
defining the volume of the pycnometer. 

The pycnometers were calibrated with freshly boiled, 
double-distilled water at 60°F and 100°F, using the values 
of the specific volume of water given in the International 


PREPARATION OF MIXTURES 


The cyclohexane and n-heptane were purchased from the 
Phillips Petroleam Company and were guaranteed by the 
manufacturer to have a purity of 99 mole % or better ; 
the benzene was Baker and Adamson’s “ reagent” grade. 


Binary mixtures of known composition were prepared 
by adding a quantity of each of the component liquids to 
a glass-stoppered bottle and weighing the bottle and con- 
tents after each addition of liquid. With the knowledge 
of the weight of the empty bottle, the weight of each of 
the liquids was then calculated and the composition 
expressed in terms of the weight fractions of the individual 


Critical Tables for these temperatures. Since the specific 
volume measurements and the calibrations were made at 
™ the same temperatures, no correction for the expansion 
of glass was necessary. 

1956/57 The weights of the empty and filled pycnometers as 
well as the weights of the pure compounds, which were 
mixed to make solutions of known composition, were 
determined with an analytical balance of 2kg capacity 
with a sensitivity of 02mg. The weights for use with 
the balance were checked against standard weights 
calibrated at the Bureau of Standards. 

The specific volumes of the pure compounds and their 
binary mixtures were determined at the two temperatures 
using different samples. In the study of the ternary 
mixtures, the specific volumes were determined with one 
loading of the pycnometer by using two thermostated 
water baths. The procedure was to immerse the pycno- 
meter first in the bath at 60°F until equilibrium was 
attained, then bring the liquid level to the reference 
mark, remove it and determine its weight ; next, transfer 
it to the bath at 100°F and, after adjustment of the liquid 
level, remove it and redetermine its weight. The water 
baths were vigorously stirred and their temperatures 
maintained to within + 0-02°F of the desired temperature 

The three liquids were used without additional purifica- 
tion. 
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components. It may be noted here that the weights of 
the components in the vapour phase were very small and 
the effect of neglecting them in the calculation of the 
weight fraction was negligible. 

Ternary mixtures were made up in a similar manner 
by blending cyclohexane in various amounts with each of 
five binary mixtures of n-heptane and benzene of constant 
composition, prepared as stock solutions and designated 
by the symbols BH1 to BH5. Compositions were expressed 
in terms of the weight fractions of cyclohexane and stock 
solution. 

Samples of these mixtures were transferred to the 
pyenometer by replacing the stopper of the stock bottle 
with a pressure siphon, similar to that used on a laboratory 
wash bottle. The end of the delivery tube of the siphon 
was drawn to a capillary so that it could be inserted 
through the capillary section of the pycnometer and the 
transfer made by applying air pressure over the liquid in 
the stock bottle. Changes in composition due to evapora- 
tion of the more volatile component during the transfer, 
were thus reduced to a minimum. 


THEORETICAL CONSIDERATIONS 


In the treatment of nonideal solutions, HiLpeE- 
BRAND [3] has defined a “ regular” solution as 
one involving no entropy change when a small 
amount of one of its components is transferred 
to it from an ideal solution of the same composi- 
tion. Such a solution is free of orienting and 
chemical effects so that the distribution and 
orientation of the molecules are random, just as 
in the ideal solution. Petroleum hydrocarbon 
mixtures may be considered as examples of 
regular solutions, since the compounds com- 
posing the mixtures are nonpolar and exhibit 
little, if any, solvation, association or chemical 
effects. 

For a ternary regular solution at constant 
temperature, H1LDEBRAND [4] gives the following 
equation for the change in volume on mixing of 
the components 1, 2, and 3: 


A V™ — (1) 
AE” V°(A,, + A, Ay (2) 
Av™ = V™— ve (3) 
B° = Bid, + Babs + Babs (4) 


4, = $3 = (6) 


and w,, w, and ws, V,, V,, and Vy, B,, By, and B, 
are weight fractions, specific volumes and iso- 
thermal compressibilities of the components 1, 2, 
and 3; A,, 4, 4, 5, and A, , are interaction con- 
stants for the respective binary pairs and V™ is 
the specific volume of the mixture. The per cent 


volume change then, is : 
1004 V™ 


4.2.3 = 100 (¢,8, + $28, + 585) 


(Ay, + Aja + (7) 


Use of this relation requires a knowledge of the 
compressibilities of the components as well as the 
interaction constants. However, by rearrange- 
ment, the equation can be put in such a form 
that the constants of the resulting equation can 
be determined in a simple manner from volume 
change data alone. 

Equation (7) in expanded form is : 


= 100 A, + + 
+ 100 A, (¢,8; + $28, + $385) 
+ 100 A, + $28. + $38). (8) 
Since ¢, + ¢, + ¢, = 1, it can be shown that, 
100 A 


100 A, (¢,8, + ¢,8, + $383) 
= bibs [100 Ay, 4) 
(10) 
+ (2 | 
The second and third terms of equation (8) can 
be similarly rewritten. Adding the three terms 
together gives : 
= [41,2 + 1,2 — 
+ [a, 5 + b, 
+ [42,5 + — $s)] 


2.3 
where, in terms of any components i, j,and k 


+ 100 *s) (¢, — 


(11) 
a; ,=1004,, (12) 


; = 100 A, ; (18) 
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2 
C123 = +S (15) 


For a binary system, on putting ¢, = 0, equation 
(11) reduces to: 


Aj; = $id; [a; ; b; (¢; — ¢;) | (16) 


Since a,; and b,;; can be determined from 
measurements of volume change in the binary 
system composed of i amd j, it follows that 
equation (11) expresses the percentage volume 
change on mixing three components in terms of 
the constants which can be derived from a study 
of the three binary systems with only the constant 
C).9.3 to be obtained from the data on ternary 
mixtures. 

Since the compressibilities of the components 
are of comparable magnitude, and furthermore, 
since the terms in “b” and “c”’ involve triple 
products of volume fractions, while the terms in 
“a” involve only double products of volume 
fractions, the contributions of the “b” terms 
and, particularly, the “cc” term will be small in 
comparison to the contribution of the terms in 
“a.” If the “ec” term is dropped, then, 


+ bib + 4 — $5) | 
+ 3 + — $3) 


and, to a high degree of approximation, the 
volume change in the ternary system can be 
expressed in terms of the volume changes in the 
binary systems. 

In the investigation of the volume change on 
mixing in ternary mixtures, the volume change 
was determined by blending a mixture of constant 
composition of two of the components, say 2 and 3, 
with varying amounts of the third component. 
The measured volume change is then, the differ- 
ence between the total volume change on mixing 
1, 2, and 3 and the volume change on mixing 2 
and 3. 

Consider the mixing of ¢,, ¢,, and ¢, volumes 
of components 1, 2, and 3, where ¢, + ¢, + ¢3=1, 
so that ¢,, ¢, and ¢, are also volume fractions in 
the ternary mixture. Let AV,, be the volume 


Ai 2.3 = + by — ¢,)] 
(lla) 
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change and A, , the per cent volume change on 
mixing ¢, volumes of 2 with ¢, volumes of 3; 
AV, 43, the volume change and A, ,, the per 
cent volume change on mixing ¢, volumes of 1 
with (¢, + ¢, + AV, ,) volumes of a mixture of 
2 and 3, and AV,, 4, the volume change and 
A, og, the per cent volume change on mixing ¢, 
volumes of 1, ¢, volumes of 2 and ¢, volumes of 3. 
Then : 


A, (17) 
= 100(AV,, + AV, 43) (18) 
(¢, + ¢3) A, s 
+ (¢, + + + AV, 5) Ai 2s (19) 


= +45)(1 + 9129) + (20) 


A, 23 + (¢, + 
+ Ay +99). (21) 
Rearranging equation (21) gives : 
A, 2.3 + $3) (22) 
($y + $5) 


Neglecting the term, O23 (4, +5) since it is 


much smaller than unity : 


A, 23 2.3 Ay 3 > ¢3)- 


Equation (22), therefore, relates the per cent 
volume change in the second step of mixing with 
the per cent volume change in the first step 
(mixing 2 and 3) and in the overall process 
(mixing 1, 2, and 3). Substituting equations (11) 
and (16) into (22a) and considering the fact that 
volume fractions in equation (16) are to be based 
on the binary system, the following equation is 
obtained for A, 4, : 


23 = +b, — ¢2)] 
+ [41,5 + 1,5 (¢, — 
bobs [42.5 + bys 
2.3 — ($2 + $5) 


($2 — | 
(2 ($2 + ¢a)) 


(22a) 


(23) 


Gg, 


Equation (28) can be rearranged to give : 
Ay 28 +4) —(¢_+¢s)}] (24) 


ogg: is given by : 


The volume of the mixed components 2 and 3, 


Table 1. Volume Changes on Mixing in the Cyclohexane-n-Heptane, Cyclohexane- Benzene, and n-Heptane- 
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+ 44) (1 + (28) 


The term ous may be neglected since it is much 


where : 
4 Ka,, smaller than unity so that : 
m 148 bes = + $5: (28a) 
+ (25) Substituting equation (28a) into (24) gives : 
Kb, Ay = [41,28 + 51,28 (¢; (248) 
1.23 "= i+K + 1+ K The resulting equation, therefore, is of the same 
Kb K form as equation (16) for a two-component system. 
+ — 33 — S61, 2.8 (26) By the same procedure the equation can be shown 
to apply to multicomponent systems of any com- 
lexity. 
4 P 
ResvuLts 


For the computation of the volume change on 
mixing in the cyclohexane-n-heptane-benzene 


Benzene Systems at 60° F 


P= Cole 


j = 


— (Aa je 
System a (Ay je | de — (Ai je | (Ai je 
| =x 100 
| 
i = 0-0 145425 | — | — 
0-18134 1-42402 0-129 + 0-005 +38 
j = n-CrHy, 035113 1-39528 0-224 0-000 0-0 
0-44672 137881 | 0-259 — 0-006 
0-63176 134634 | 0-284 — 0-009 — 32 
0-83535 130927 0-211 + 0-005 +24 
1-27739 


1-138220 


1:15618 0-333 + 0-007 + 21 
1-17952 0-547 — 0-012 — 22 
110444 0-647 — 0-001 — 02 
1-21317 0-684 — 0-009 —13 
124038 0-599 + 0-008 +13 


127739 


1:13220 


1-20080 0-377 + 0-009 + 2-4 
1-26726 0-508 — 0-012 —24 
1-29639 0-521 — 0-006 —12 
130434 0-518 — 0-005 —10 
1:31949 0-506 — 0-002 — 04 
1:89580 0-306 + 0-005 +16 


145425 


1956 


VOL 

0-0 | | = | 
0-13886 
0-28188 
0-87588 
0-50085 
0-69422 
i= 0-0 | _ | --- | 
0-19903 
0-48893 
0-51361 
056086 
0-80535 
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system, the specific volumes of the 3 pure com- 
ponents, 5 mixtures of each of the 3 binary 
systems and 25 mixtures of the ternary system 
were determined at 60°F and 100°F. These data 
are given in Tables 1, 2, 3 and 4. 

All of the data on the binary mixtures and the 
data on the ternary mixtures prepared by blending 
binary mixtures of constant composition of 
n-heptane and benzene with cyclohexane, were 
treated by application of the equations (16) and 
(24a) respectively. For this purpose both equa- 
tions were rewritten in terms of volume change 


= 4; + (¢; $;) (29) 
where i is a pure component and j may be either 
a pure component or a binary mixture of constant 


composition. Equation (29) may be put in the 
form : 
100 AV™ 


m, ; = (a, ; — 5; ;) (31) 


= ; nj (30) 


where 


and 
= 2 b; (32) 


By equation (30) it is evident that a plot of 
the volume change coefficient vs. ¢; will result 
in a straight line with an intercept of m;; and 
a slope of n; ; from which the values of a; ; and 
b, ; can be obtained. Figs. 1, 2, 3, and 4 show 
the experimental data plotted in this manner. 
Best straight lines were drawn through the points 
for a given system in order to evaluate the con- 
stants a; ; and 6; ; which are listed in Table 5. 


Using these values and equations (16) and (24a) 


Table 2. Volume Changes on Mixing in the Cyclohexane-n-Heptane, Cyclohexane- Benzene and n-Heptane- 
Benzene Systems at 100° F 


( Mi je ( Ai je 


x 100 


: 100 AV™ 
coefficient, as: 
System (Ay je (Aj de — (Ai, (Aj, je 
0-17838 146442 0-073 — 0-004 — 55 
j = n-CrHy, 0-36396 1-43200 0-157 + 0-007 +45 
0-47574 1-41206 0-180 — 0-002 -—1}1 
0-60884 1-38811 0-195 — 0-004 —20 
0-78944 1-35495 0-166 0-000 0-0 
1-0 131458 om 
i = 0-0 1-16422 
0-23812 1-20544 0-451 + 0-001 +02 
j = CoH, 0-41788 1-23455 0-612 + 0-001 + 02 
0-51349 1-24932 0-637 + 0-006 +09 
0-68186 1-27373 0-552 + 0-002 +04 
0-85576 1-29691 0-311 — 0-002 +06 
1-0 1-31458 one 
0-18079 1-22789 0-306 — 0-002 — 06 
j = CoH, 0-35679 1-28827 0-451 — 0-008 — 06 
0-51125 134009 0-480 + 0-008 +17 
0-66142 1-38912 0-410 + 0-005 +12 
0-80390 1-43446 0-265 — 0-008 — 30 
15 
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VOWME FRACTION N-HEPTANE 


VOLUME CHANGE COEFFICIENT 


0204 06 08 
VOLUME FRACTION CYCLOHEXANE 6. 
Fic. 1. Effect of composition on volume change coefficient Fig. 2. Effect of composition on volume change coefficient 

for binary mixtures at 60°F. for binary mixtures at 100°F. 


VOLUME FRACTION CYCLOHEXANE 9, VOLUME FRACTION CYCLOHEXANE @, 


Fic, 3. Effect of composition on volume change coeflicient Fic. 4. Effect of composition on volume change coefficient 
for ternary mixtures at 60°F. for ternary mixtures at 100 F. 
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the curves of per cent volume increase vs. weight 
fraction, shown in Figs. 5, 6, 7, and 8 were 
calculated. The experimental values, represented 
by the circles, have been added for the sake of 
comparison. Values of a, ,, and b, ,, for the 


mixing of cyclohexane with stock solutions BH1 
to BH5 were calculated by means of equations 


(25) and (26) using the values of the constants 
listed in Table 5. These calculated values are 
compared with the experimental values obtained 
by equation (24a) in Table 6. 

In order to obtain the volume changes on mix- 
ing the pure components cyclohexane, n-heptane 
and benzene, it was necessary to know the volume 


Table 3. Volume Changes on Mixing Cycloherane with Miztures of n-Heptane and Benzene at 60°F 


y™M 


(Aj, je 
x 100 


(Aide (M4 se — 


i = 
j = BHe 


K = 0-4130 


i = 0-0 139438 t 
0-21562 1-37147 0-148 + 0-002 +14 
j= BHI 0-33153 1-35875 0-201 — 0-004 —20 
VOL. 0-45450 1-84513 0-249 + 0-002 +08 
6 K = 0-1939 0-65202 132227 0-252 — 0-002 ~08 
956/57 0-78738 1-30597 0-206 + 0-004 +19 
0-14614 1-34004 0-109 — 0-005 ~ 46 
0-30028 1-33059 0-212 + 0-004 +19 
0-46105 1-32003 0-269 + 0-001 + 0-4 
0-65002 1-30668 0-265 0-004 
| 0-79518 1-29577 0-210 + 0-004 +19 
10 1-27868 os 
i = 0-0 1-30023 
j = BH3 0-19705 1-29863 0-204 + 0-005 +24 
0-31687 1-29697 0-276 0-008 =~ 34 
K = 0-7768 0-47829 1-29423 0-334 0-000 0-0 
| 0-65214 1-29021 0-314 — 0-001 —03 
0-78416 1-28644 0-242 0-000 0-0 
10 1-27868 on 
i = 0-0 1-25152 
j = BHs 0-20849 126032 0-249 + 0-008 +33 
0-33717 1-26493 0-337 0-004 
K = 1-402 0-50135 1-27011 0-393 0-008 
0-67132 1-27450 0-374 + 0-008 +08 
0-80673 1-27686 0-269 — 0-008 — 13 
10 1-27868 ~ 
i = 0-0 1-20167 
j = BHS 0-20364 1-22091 0-292 0-002 0-7 
0-31676 1-28102 0-404 + 0-008 +07 
K = 3-0899 0-51501 1-24781 0-482 0-000 0-0 
0-64664 1-25710 0-450 ~ 0-001 ~02 
0-79876 1-26728 0-324 — 0-001 ~08 
17 
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changes occurring in the preparation of the stock 
solutions BH1, BH2, BH3, BH4, and BH5. From 
the specific volumes of the stock solution and the 
equation for the volume changes in the benzene- 
heptane system, established earlier, the exact 
compositions of these solutions were calculated. 
Hence, the per cent volume changes that occurred 


Table 4. 


in the preparation of the stock solutions as well 
as their ideal volumes could be calculated. Equa- 
tion (21) was then applied to calculate the per 
cent volume change, A, ,,, from the experi- 
mental data. The values are shown along with 
the compositions of the ternary mixtures in 
Tables 7 and 8. 


Volume Changes on Mixing Cycloherane with Mixtures of n-Heptane and Benzene at 100° F 


System Wi 


(Ai je 


(Si. — (Ai, je 
(Aj, je 
100 


(Ai, — (Ai, je 


00 1.43349 
0-21562 140968 
033153 139645 
045450 138225 
0-65202 1:35859 
0-78738 134186 
10 131385 


138637 
137734 
136732 
1-35631 
134259 
133131 
1:31385 


133633 
133456 
133271 
132975 
132562 
132186 
131385 


128616 
120519 
129981 
130515 


020849 
033717 
0-50135 
0-67 132 130952 
080673 131197 
10 131385 


123473 
125464 
126499 
128170 
129167 
130216 
131385 


0-004 
0-005 
0-004 
0-007 
0-005 


0-009 
0-007 
0-003 
+ ool 
+ 0-002 


j = BHI +29 
0-188 - — 26 VOL 
K = 0-1939 0-227 - —17 6 
0-229 30 195€ 
0-192 + +27 
= BH2 0-14614 + 0-002 +18 
030028 0-199 0-003 —15 
K = 0-4130 046105 0-249 — 0-007 —27 
0-65002 0-250 0-003 —12 
0-79518 0-196 + 0-004 +21 
10 
i = C,H, 0-0 
j = BH3 019705 0-200 0-009 +47 
031687 0-264 0-004 — 15 
K = 0-7768 047829 — 0-006 —19 
0-65214 0-299 0-005 —16 
O-78416 0-239 + 0-006 +26 
10 
j= Bs 0-253 + 37 
0-334 —21 
K = 1-102 0-394 
| 0-366 +31 
0-266 + 
i= | 00 | 
j = BHS 020364 0-304 + 0-001 + 
031676 0-413 + 0-002 + O05 
K = 30899 0-51501 0-487 0-000 0-0 
064664 0-449 — 0-003 
T9876 0-326 | + 0-003 +09 
| | 
18 
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Fic. 5. Effect of composition on volume change in binary 
mixtures at 60°F. 


A, was also calculated by equation (11) 
and by the approximate equation (11a) using the 
constants a, and c,,, given in Table 5. 
2,9 Was determined by the method of least 
squares from the difference between the values 
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Fia. 7. Effect of composition on volume change in ternary 
mixtures at 60°F. 
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Fic. 6. Effect of composition on volume change in binary 
mixtures at 100°F. 


of A,,, calculated by equation (21) and 
equation (lla). The resulting values of A, , , 
are given in Tables 7 and 8. By the use of a 
digital computing machine to solve equation (11) 
a series of points of constant per cent volume 
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Fic. 8. Effect of composition on volume change in ternary 
mixtures at 100°F. 
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Table 5. Experimental Values of Constants of Equations (16) and (11) at 60°F and 100°F 


System 60°F 100° F 
i j j bi, j | j | 
| | 
1-150 0-350 0-770 0-290 
| CoHy | 2-760 0-380 2510 0-360 
2-110 0-210 1-900 0-200 
BHI | 1-035 0-215 0-965 0-185 
BH2 | 1-105 | 0-215 1-050 0-160 
BH3 1-345 0-125 1-295 0-125 
BH4 1-600 0-270 1-585 0-215 
BHS5 1-905 0-285 1-040 


Values of C, 9 5 of Equation (11) are: — 1-00 at 60°F and + 1-40 at 100°F 


Table 6. Comparison of Experimental and Calculated Interaction Constants for Mixing Cyclohexane 
with Mixtures of n-Heptane and Benzene at 60° F and 100° F 


i = 
j = BHI 


K = 0-1939 


= 
j = BH2 


K = 0-4130 


i = 
j = 


K = 0-7768 


i = 
j = BH 


= 
j = BHS 
K = 8-0899 


1-035 
1-013 


1-345 
1-176 


1-600 
1-457 


1-905 
1-895 


100° F 
0-215 0-965 0-185 
0-447 0-848 0-229 
O215 1-050 0-160 
0-493 0-985 0-194 
0-125 1-295 0-125 
0-515 1-202 0-176 
0-270 1-585 0-215 
0-510 1-500 0-182 
0-285 1-940 0-240 
0-474 1-873 0-223 


increase were computed for the ternary system. 
These data are shown graphically (by the curves) 


in Figs. 9 and 10. 


Precision oF RESULTS 
The experimental procedure was kept uniform as 


far as possible so that the accuracy of the volume 
change data would be independent of any 
systematic errors in the individual specific 
volume measurements. The limits of precision 
of the various operations involved and the 
corresponding uncertainty in the specific volume 
measurements were estimated to be as follows : 


vol 
1956 
System 
9}, 23 
K = 1-4402 
20 
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Fic. 9. Lines of constant per cent volume increase on mixing in the cyclohexane-n-heptane-benzene system at 


60°F. 


Fic. 10. Lines of constant per cent volume increase on mixing in the eyclohexane-n-heptane-benzene system at 
100°F. 
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Operation 

Weighing 

Temperature 

Adjustment of liquid level 


in pycnometer 0-2 mm 


Check determinations of the specific volume 
were found to be reproducible to within 8 parts 
in 100,000. 


Discussion 


The relations between the composition and the 
per cent volume change on mixing in the cyclo- 
hexane-n-heptane-benzene system at 60°F and 
100°F, are summarized in Figs. 9 and 10. It is 
to be noted that throughout the entire range of 


composition an expansion in volume takes place 
on blending, with a maximum volume change of 
about 0-7% occurring somewhere in the region 
of a 45-10-45 cyclohexane-n-heptane-benzene 
mixture. Of the binary systems at 60°F, the 
cyclohexane-benzene system shows the greatest 
expansion with a maximum of about 0-7% and 
the cyclohexane-n-heptane the least with a 
maximum of 03% while the heptane-benzene 
system is intermediate with a maximum of about 
0-54%. Increasing the temperature to 100°F 
reduces the values somewhat. 

That the sign and the magnitude of these 
volume changes are to scme extent associated with 
differences in the chemical nature of the com- 
ponents is substantiated by mixing data on these 
same compounds with other naphthenic, paraffinic 


Table 7. Volume Changes on Mixing in the Cyclohexane-n-Heptane-Benzene System at 60°F 


2, ae 


2, ale — (4i, 2, ade 
2, ale 
x 100 


(44, 2, ale — (41, 2, ade 


Eqn. (21) 


Eqn. (lla) | Eqn.(11) | Eqn.(1la) Eqn. (11) 


— 0-006 
— 0-019 
— 0-022 
— 0-033 
— 0-027 
— 0-013 
— 0-015 
— 0-027 
— 0-038 
— 0-025 


+0011 
+ 0-001 
— 0-003 
— 0-021 
— 0-021 
+ 0-009 
+0015 
+ 0-001 
— 0-021 
— 0017 
+ 0-041 
+ 0-082 
+ 0-019 
— 0-006 
— 0-016 
+ 0-026 
+ 0-014 
+ 0-002 
+ 0-001 
— 0-009 
— 0-008 
— 0-004 


— 15 
— 46 
— 52 
— 90 
— 26 


+28 
+ 0-2 
— 07 
— 57 
— 7-7 
+18 
+28 
+02 


— 0-004 


— 0-026 
— 0-025 
— 0-006 
— 0-022 
— 0-028 


Uncertainty in 
Specific Volume 
Precision parts per 100,000 
05 mg 12 
0-02 °F 20 
08 
Vol 
6 
| 195¢ 
weight fraction 
| 
CoH, 
| 
0-2156 0-6280 0-1564 0-393 | | 
0-83815 0-5852 0-412 
0-4545 0-4367 0-1088 0-423 
0-6520 0-2786 0-0694 0-365 
0-7874 0-1702 0-0424 0-273 
0-1461 0-5579 0-2960 0-501 
0-30038 0-4571 0-2426 0-536 —28 
0-4611 0-3521 0-1868 0-521 — 52 
0-6500 0-2287 0-1213 0-430 —88 —49 
0-7952 0-1838 0-0710 0-307 —81 —55 
0-1970 0-4018 0-4012 0-629 + 0-010 +16 +65 
0-3169 0-3418 0-3413 0-639 FP — 06 +50 
0-4783 0-2611 0-2606 0-612 — 0-013 -—21 +31 
0-6521 0-1741 0-1738 0-499 — 52 —12 
0-7842 0-1080 0-1078 0-357 —70 —45 
0-2085 0-2777 0-5188 0-646 —09 +40 
0-83872 0-2825 0-4808 0-669 — 33 +21 
0-5014 0-1748 0-8238 0-642 —44 +03 
0-67138 0-1158 0-2134 0-537 -— 0-016 —380 +02 
0-8067 0-0678 0-1255 0-365 — 0-016 —-44 — 255 
0-2086 0-1602 0-6362 0-585 — 0-027 — 46 —05 
0-8168 0-1875 0-5457 0-654 — 0-081 —-47 — 06 
0-5150 0-0976 0-8874 0-657 — 0-088 — 0017 — 58 —26 
0-6466 0-0711 0-2823 0-576 — 0-085 — 0-021 — 61 — 36 
0-7988 0-0405 0-1607 0-395 — 0-024 — 0-019 — 60 —48 
22 
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and aromatic compounds reported in the litera- 
ture. In a study of the effect of differences in 
chemical nature and physical properties of 
hydrocarbon compounds on their volume changes 
on mixing, Iorre [5] found that, in general, 
when aromatic hydrocarbons are mixed with 
nonaromatic hydrocarbons an expansion occurs 
and that the expansion is greater with naphthenic 
than with paraffinic compounds. On the other 
hand, when paraffinic compounds are mixed with 
one another, volume contractions occur. Kurtz 
and Sankin [6] and ScuiessLer and co-workers 
{12] found that only small volume changes, 
frequently contractions, occur on mixing non- 
aromatic compounds whereas large volume 
expansions are obtained in mixtures containing 
aromatics. These observations are in agreement 
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with what has been found in this laboratory in 
the blending of petroleum distillates, namely 
that straight run Mid-Continent distillates (rich 
in paraffinic compounds) show contractions when 
blended with one another and expansions when 
blended with hydroformate distillates (rich in 
nonparaffinic compounds) [8]. 

In general, satisfactory agreement is obtained 
between the experimentally derived volume 
changes and the values calculated from the 
equations for regular solutions. Inspection of 
Tables 1, 2, 3, and 4, and Figs. 1 through 8 shows 
that the average deviation of the experimental 
per cent volume change for the binary systems 
from that calculated by equation (16) amounts 
to + 0-004 units or 1-5%, which is of the same 
order of magnitude as the estimated experimental 


Table 8. Volume Changes on Mixing in the Cyclohexane-n-Heptane- Benzene System at 100° F 


(A), 2, de. — (4, 2, se 


Composition (41, 2, ae 
weight fraction 
Eqn. (21) 
0-3315 0-5352 0-1333 0-378 
04545 0-4867 0-1088 0-384 
0-6520 0-2786 0-0694 0-331 
0-7874 0-1702 0-0424 0-255 
0-1461 0-5579 02960 0-467 
0-3003 6-4571 02426 0-491 
0-4611 0-3521 0-1868 0-476 
0-6500 0-2287 01213 0-399 
0-7952 0-1338 0-0710 0-283 
0-1970 0-4018 0-4012 0-582 
0-3169 0-3418 0-3413 0-589 
0-4783 0-2611 0-2606 0-565 
0-6521 0-1741 0-1738 0-466 
0-7842 0-1080 0-1078 0-343 
0-2085 0-2777 05138 0-609 
0-3372 0-2325 0-4303 0-631 
0-5014 0-1748 0-3238 0-616 
0-6713 0-1153 0-2134 0-512 
0-8067 0-0678 0-1255 0-351 
0-2036 0-1602 0-6362 0-565 
0-3168 0-1375 05457 0-636 
0-5150 0-0976 0-3874 0-644 
0-6466 0-0711 0-2823 0-562 


(Ay, 2, ale — (41, 2, ale (A; 2, ve 
x 100 
Eqn. (11a) | Eqn.(11) | Eqn.(lla) | Eqn. (11) 

+ 0-036 +0012 + 10-0 +33 
+ 0-035 + 0-007 +92 +19 
+ 0-037 +0014 +96 + 3-6 
+ 0-025 + 0-009 +75 + 2-7 
+ 0-024 + 0-016 +94 +63 
+ 0-024 — 0-006 +51 —13 
+ 0-081 — 0-012 +63 —2-4 
+ 0-027 — 0-013 +57 —2-7 
+ 0-021 — 0-003 +53 —08 
+ 0-016 + 0-005 +57 +18 
+ 0-046 + 0-002 +79 +03 
+ 0-039 — 0-012 + 6-6 —20 
+ 0-037 — 0-008 +65 —14 
+ 0-025 — 0-003 + 54 — 0-6 
+ 0-023 + 0-010 +67 +2-9 
+ 0-039 — 0-006 + 64 —10 
+ 0-030 — 0-020 +47 — 3-2 
+ 0-036 — 0-006 +58 —10 
+ 0-035 +0011 + 68 +21 
+ 0-023 + 0-013 + 66 
+ 0-023 — 0011 +41 —19 
+ 0-026 — 0-012 +41 -—19 
+ 0-030 0-000 +46 0-0 
+ 0-022 + 0-002 +39 + 0-4 
+0019 +49 +28 


error. Table 6 shows that reasonably good 
agreement is obtained between experimental and 
calculated values of the constants a, ,, and b, 9 
for mixing cyclohexane with constant composition 
mixtures of n-heptane and benzene. Tables 7 
and 8 show that the approximate equation (Ila), 
based on binary data alone, predicts the ternary 
volume changes to within about + 5% of the 
experimental value. The use of a ternary inter- 
action constant, ¢c, , 5, brings the agreement to 
within + 2-5%. 

The per cent volume changes on mixing in the 
cyclohexane-benzene and n-heptane-benzene 
systems are in reasonably good agreement with 
the data reported on these systems in the litera- 
ture. Woop and Austin [12] report per cent 
volume changes in a 50-50 mixture by weight of 
cyclohexane and benzene of 0-661 and 0-657 at 
60°F and 100°F respectively, whereas 0-693 and 
0-631 were found in this research. Brown and 
Ewa.p [1] report a per cent volume change in 
a 50-50 mixture of n-heptane and benzene at 25°C 
(77°F) of 0-531 which may be compared to 0-526 
at 60°F and 0-474 at 100°F reported here. Kurtz 
and Sankin [6] report 0-27% volume change in 
n-heptane-cyclohexane mixture containing 53-2% 
n-heptane at an unspecified temperature while 
0-287% volume change was obtained in the 
present study for the same composition at 60°F. 


CONCLUSIONS 


It ha: been shown that the volume expansion 
that occurs when cyclohexane, n-heptane and 
benzene are mixed at room temperature, can be 
appreciable (0-7%) in certain concentration 
regions. It was shown also that the volume 
change on mixing benzene, an aromatic, with 
either cyclohexane, a naphthene, or n-heptane, 
a paraffin, is much larger than that on mixing 
cyclohexane with n-heptane. It is concluded that 
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The separation of multicomponent mixtures of common gases by 
thermal diffusion 
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Abstract—A brief resumé is given of theoretical and experimental work on the separation of 
binary and multicomponent mixtures of isotopes. The use of this information for the design 
of vertical hot-wire types of thermal diffusion columns is discussed. It is pointed out that no 
comparable experimental work exists for the separation of multicomponent common gases. The 
existing theory based on isotope separation is inadequate when dealing with common gases. 


Experimental results are quoted for total reflux and continuous operation of a pilot-plant 
thermal diffusion column consisting of concentric tubes, and working on common gas mixtures. 
It is demonstrated from the results, that the separation of CO from Hg may be enhanced by 
the introduction of CO, which produces a maximum separation of CO from Hy. The possible 
reasons for the effect are discussed. 


In conclusion, the absence of theory and the necessary complexity of such theory, when 
developed, indicates the necessity for a semi-empirical approach to the design of thermal diffusion 
columns for separation of multicomponent mixtures, based on pilot-plant information. 


Résumé— Les auteurs donnent un rapide résumé d'un historique sur la séparation théorique et 
expérimentale des mélanges d'isotopes binaires ou & plusieurs constituants. Ils discutent usage 
de ces données pour létude des colonnes de diffusion thermique 4 fil chauffant vertical. Insistant 
sur le fait qu'il n’existe pas de guide pour la séparation de gaz usuels a plusieurs constituants, ils 
appliquent & ces gaz la théorie de la séparation des mélanges binaires d'isotopes. 

Les auteurs citent les résultats expérimentaux obtenus avec une colonne pilote de diffusion 
thermique composée de tubes concentriques marchant en continu et a reflux total pour des 
mélanges de gaz usuels. Les résultats montrent que la séparation entre CO et H, peut étre 
facilitée par l'introduction de CO,. Les causes possibles de cet effet sont discutées. 


L’absence de théorie et la complexité obligatoire d'une telle théorie, quand elle est développée, 
montrent la nécessité d'opérer semi-empiriquement pour le calcul d'une colonne de diffusion 
thermique pour la séparation des mélanges multiples en se basant sur les données d' installations 
pilotes. 


INTRODUCTION 
Tue difficulty of solution of problems on the 
design of thermal diffusion apparatus from the 
separation of binary and/or multicomponent 
mixtures of common gases can be best understood 


The theory of separation of binary miztures of 
isotopes 
The classical theory is due to Cuapman [1] and 


Enskoc [2]. Their treatment of the subject is 
extremely complex, and various other theories 


by a review of the information available on isotope 
separation. The results of an experimental 
investigation on the separation of ternary mix- 
tures of common gases in thermal diffusion pilot 
plant will then be discussed in the light of the 
review. 


with a more elementary approach have been 
proposed by [3], Franxer [4], and 
Forrnu [5]. The theory of Franke and Jones, 
developed by Fixrn, has been extended by 
Cacctapuoti [6], Furry [7], and Furry and 
Jones [8]. 
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A simple treatment has been given by WHALLEY 
and WINTER [9], essentially based on the work 
of Firtu. Their solutions are not exact, but 
have a value for first analysis. Theoretical values 
of the thermal diffusion coefficients k; and « were 
obtained by Winter for mixtures H, - He, 
H,-N,, H,- Ne, N,-CO,, and they showed 
that there was a marked temperature dependence 
of these coefficients not allowed for in the theory. 
Further, it was evident that the agreement 
between values of the coefficients calculated by 
exact and elementary theories depended on the 
binary mixture being considered, and in some 
cases lack of agreement was considerable. 


The theory of separation of multicomponent 
mixtures of isotopes 

The general theory for such separations has been 
given by and Uex.ine [10], but their 
equations are very complex and difficult of 
application. and Winter [9] have 
extended their elementary approach used for 
binary mixtures, to the case of multicomponent 
systems, and in particular ternary mixtures of 
isotopes. They considered the effect of He on 
the separation of Ne isotopes, and also, the 
effect of Ne on the He isotopes. Their calculated 
values are not exact, but they indicate that the 
magnitude of the effect of the third gas is not 
large. 

WINTER investigated the results of Cu1pMAN 
and Dastvr [11], and first published experimental 
work on this subject. They introduced an inert 
gas into a steam-hydrogen mixture, and studied 
its effect on the apparent equilibrium constant 
of the system. 

H, (g) + O (in Fe) = H,0 (g) 
While Winter accepted a qualitative agreement 
with his equation, he could not draw any 
quantitative conclusions from the work. 


The separation of multicomponent mixtures in 
vertical hot-wire type columns 

The first published work is that of E_per [12]. 
It is a brief report of an experimental investigation 
on the separation of the neon isotopes with mass 
20, 21, and 22 respectively. The apparatus con- 
sisted of a 2-stage 6m long. glass column of 


The separation of multicomponent mixtures of common gases by thermal diffusion 


hot-wire type, containing 38 litres of gas. Total 
reflux operation only was carried out, the con- 
centrations at the end of the equilibrium period 
being recorded. The reported experimental 
operating data are inadequate to form an opinion 
as to the effect of any one component on the 
separation of the other two. 

This work was extended by ELper and 
McDonatp [13], the separation of the neon 
isotopes again being considered. The apparatus 
consisted of a 3m glass tube of hot-wire type. 
The operating conditions were not well defined 
in the report. Compositions along the column 
were obtained at selected points by means of a 
mass spectrometer. On a plot of concentration 
versus height for each of the isotopes, it was 
found that the distribution of neon 21 along the 
column was not a smooth monotonic curve as for 
neon 20 and 22. This they attributed to a partial 
differential separation of neon 21 and 22, which 
is superimposed on the predominating separation 
of both of the heavier isotopes from 20, The result 
does not give a positive indication of the effect 
of the introduction of a third component on the 
separation of a binary mixture. 

In 1953, Ciusitus and Scnumacner [14] 
examined the problem of the separation of poly- 
mixtures of isotopes in CLusius-DickEL type 
columns. A comprehensive theoretical treatment 
based on CHapMAN [1], JENSEN [15] and HeLLtunpD 
[10], is given for binary mixtures, and poly- 
mixtures. An assumption is made that the 
enhancement of the separation produced by a 
vertical column is independent of the number of 
components present. This may be true as long 
as the density remains virtually independent of 
the concentration along the column, and would 
apply reasonably to isotope mixtures. It is 
assumed but not stated, that k,and « the thermal 
diffusion coefficients, are independent of con- 
centration and temperature. SCHUMACHER claims 
that in mixtures of isotopes, the concentration 
of each component may be considered on its 
instantaneous value, and is not determined by 
other components. The solutions obtained are 
not exact, even for three component mixtures. A 
three component system is considered where a 
third component acts in such a way as to improve 
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the separation of the other two components over 
that separation which otherwise would be 
obtained. The choice of the third gas is based on 
a consideration of symmetry of the molecules, 
there being a maximum effect when the “* Hilse- 
gasse’ has a mass corresponding to the mean 
mass of the other two components. This gas 
would concentrate at some definite height in a 
column, and so would act as a buffer zone between 
the light and heavy ends, being easily separated 
from either. Krypton and xenon were separated 
in this way using silicon tetrafluoride as the third 
component, Extremely small withdrawal rates 
were used. The isotopes of argon were separated 
in continuous separation, but again long periods 
were allowed for equilibrium separation to be 
established. The theory was used to predict the 
size of apparatus required, but for krypton and 
xenon this could only be very approximate, and 
indeed it is doubtful whether the theory could 
be applied generally. 

In conclusion, it is apparent that all existing 
theoretical and experimental work applies to 
either isotopes or cases where extremely small 
withdrawal rates are maintained. All the appara- 
tus is of hot-wire type, with an extremely small 
tube. Such methods would be of little use for an 
industrial separation of common gases from 
multicomponent mixtures. 

The difficulties inherent in the design of 
thermal diffusion columns for the separation of 
binary mixtures of common gases, have already 
been examined by the authors [16]. In 1953, 
work was carried out on a single thermal diffusion 
column of pilot plant dimensions. The experi- 
mental study was of a preliminary nature only, 
and was a logical development of the work carried 


out in this apparatus on binary mixtures. The 
work is discussed below. 


An EXPERIMENTAL StTupy OF THE 

Separation BY TuermMat DIFFUSION 

or A Turee ComMPONENT MIXTURE OF 
Common GASES 


The gases chosen were H,, CO,, and CO. The 
first experiments were to determine how the 
introduction of CO, affected the separation of 
H, from CO, and secondly, how CO affected the 
separation of H, from CQ,. 

The apparatus used and the operating tech- 
niques employed have already been described 
by the authors [16]. For the present experiments 
a 6ft. column made in steel was used, com- 
prising concentric cylinders. The dimensions of 
this column are given below. (This is column 2, 
described in the authors’ paper ref. [16], where 
the inner tube is of stainless steel. This is not 
good practice due to enhancement of temperative 
asymmetry. However for column 2, it is felt 
that the effect is not so pronounced as to affect 
appreciably the results.) 

The column was operated both under “ total 
refiux ", and “ continuous” conditions. In the 
present case, the total reflux operation refers to 
the column closed at the lower end, but with 
circulation of gas in a closed circuit around the 
top of the column, a large gas holder being used 
to ensure a substantially constant top con- 
centration. The continuous operation on the 
single column, refers to top circulation of gas, 
but with a continuous withdrawal of gas from 
the lower end of the column. Further, for common 
gases, this withdrawal rate is a substantial one 
as will be seen from the experimental figures. 


EXPERIMENTAL RESULTS 


Table 1. 


Column dimensions, column 2. 


Inner diameter of outer column 
Outer diameter of inner column 


Heated length 
Annulus width 


6 ft. 


Mild steel outer tube. Stainless steel inner tube. 
Electrically heated element within the inner tube. 
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of CO, into H,/CO mixtures. 


Table 2. Total reflux column 2, T, = 270°C, AT = 250°C, 2w = 180m, Introduction of CO, 
into H,/CO miztures. 

CO, | co | H, | H,/CO | P | Ac, 
Initial gas 6-2 25-6 62 | 278 | 26-9 0 
Equilibrium end concentration 120 51-0 37-0 0-74 57-5 30-6 

| 
Initial gas 560 6190) || 0 
Equilibrium end concentration 2700 | 088 58-0 30-0 
Initial gas 272 562 | B47 0 
Equilibrium end concentration 38-0 28-0 35-0 | 1-20 445 23-9 
Initial gas 43-0 | 12-0 45-0 3-75 21-1 0 
Equilibrium end concentration 58-0 17-0 25-0 1-47 WS 19-4 
Initial gas oo 85-0 | 566 (150 0 
Initial gas oo | 220 | wo | 856 | 220 0 
Equilibrium end concentration 00 619 | | 16-1 
P = Percentage CO to a basis of CO, free gas mixture. 

Table 8. Continuous operation, column 2. T, = 270°C, AT = 250°C, 2w = 1:80 cm. Introduction 


Initial gas concentration 
Final gas concentration at 
o = 480 cc/min 
o = 840 cc/min 


Hy 


H,/CO 


P 


2-73 


1-24 


Initial gas concentration 

Fina] gas concentration at 
o = 860 cc/min 
= 670 cc/min 

@ = 1,000 ce/min 


8-37 


1-57 
1-77 
2-02 


Initial gas concentration 
Final gas concentration at 
o = 430 cc/min 

o = 860 cc/min 


8-47 


1-74 
2-18 


Initial gas concentration 
Final gas concentration at 
= 76 cc/min 
o = 300 cc/min 


50-6 
53-0 


12-0 


33-9 
30-7 


3-75 


2-20 
1-88 


| 


31-4 
36-3 


P = Percentage CO, to a basis of CO free gas mixture ; ¢ = withdrawal rate cc/min. 


VOL, 
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(956/57 

CO, co | | Ac 

6-2 25-6 | | | 26-9 0-0 

10-0 or | 17-7 

9-0 36-0 55-0 | 157 800 12-1 

| 

15-0 19-0 66-0 — | 23-0 0-0 

21-0 82-0 47-0 | 40-5 17-5 

19-6 29-0 51-4 | 96-1 13-1 

19-0 26-8 54-2 33-1 10-1 

27-2 | 16-2 56-2 | 22-4 0-0 

37-0 23-0 0-0 | | 36-5 145 

34-4 20-6 45-0 | 31-4 9-0 

43-0 45-0 21-1 0-0 

| | | | 

16:3 15-2 

29 
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Table 4. Total reflux, column 2, T, = 270°C, AT = 250°C, 2w = 1-80 ems. a eee 
H,/CO, mixtures, 


C0, co | H, H,/CO, 


Initial gas concentration 18-0 20-8 61-2 3-4 
Equilibrium end concentration 31-0 36-8 33-0 1-06 


Initial gas concentration 29-6 + 69-0 | 2-38 
Equilibrium end concentration 53-6 45-9 0-86 


Initial gas concentration 26-0 . 66-6 2-56 
Equilibrium end concentration 45-6 36-8 0-81 


Initial gas concentration | 22-0 78-0 | 3-54 
Equilibrium end concentration 41-2 58-8 1-43 


Q = Percentage CO, to a basis of CO free gas mixture. 


Gas-hoider 


Fic. 1. Thermal diffusion pilot plant (detailed description Ref. [15]). N — needle valves; S — sample cocks; Tr.C - 
three-way cocks. (N.B.- Only Column 2 was used for the experiments on the separation of multicomponent 
mixtures.) 


The system employed is represented in Fig. 1, 
and Fig. 2 describes the column. Gas analyses 
were obtained with a katharometer, and by 
Orsat analysis. The temperature of the inner 
tube was measured continuously, and also that 
of the cooling water. Care was taken to maintain 
a balanced operation. All experiments were 
conducted at atmospheric pressure. 

In making up the gas mixtures in the gas 
holder at the beginning of a new set of experi- 
ments, it was difficult to predict the separation, 
and hence to reproduce the H,/CO and H,/CO, 
ratios when steady-state conditions had been 
achieved. However they were obtained sufficiently 


closely to indicate any effects. Fic. 2. 6 ft. column (No. 2). (Detailed description Ref. 
Gas rates were measured by a calibrated (15}). ;, 


22-8 0-0 
48-5 25-7 
29-6 0-0 
54-0 24-4 
28-1 0-0 
| 55-4 27-3 
220 0-0 
| 412 19-2 
VOL 
6 
N ae N 
| 
s 
N 
s s g | 
I 
Heating element 
Original cooling water 
circulating —-G |i 
gos | 
| steer 
-stainiess stee! 
| 
| 
product gas 
Drain 
30 


VOL. 
6 


956/57 


rotameter on the top circulating gas, and by a 
direct bubble flowmeter on the withdrawn gas. 

In the case of total reflux, equilibrium was 
established in one or two hours. 


Discussion 
An attempt has been made to predict the results 
obtained for the total reflux condition from the 
theory as developed by CLUsius and ScuuMACHER. 
This was not successful, and it would appear that 
equations based on the separation of isotopes do 
not apply to common gases. 

The experiments on ternary 
common gases recorded in this paper. appear to 
be the first carried out in pilot plant. 

Ever carried out experimental work with a 
ternary mixture of neon isotopes. Although he 
obtained an inflection in the Ne curve, the con- 
centrations of the heavy component continued to 
increase down the column. This is also true of the 
separations carried out by Ciusivs and Scuv- 
MACHER. In the latter case. it was established 
that the third component of intermediate pro- 
perties to the pair of isotopes being separated 
tended to concentrate somewhere intermediate 
between the top and bottom of the column. 

No samples were taken along the tube in the 
present investigation, but it is doubtful if the 
third component, either CO or CO,, acts in a 
similar way in hydrogen. The molecules differ 
considerably from each other, and neither can 
be considered as being intermediate in properties 
between hydrogen and the other component. It 
would have been more accurate if the same 
initial concentration of the heavy component 
could have been used, but as will be seen from 
the H, CO and H, CO, ratios, the deviation is 
not very large, and the result obtained based on 
the enrichment ¢ is large enough to accept in 
principle. 

It will be observed from Fig. 3 for total reflux 
operation, that when CO, is introduced into a 
H, CO mixture, the enrichment of CO in H, at 
the bottom of the column, passed through a 
maximum. The initial concentration at the top 
of the column remains virtually enchanged by 
the use of a large gas holder in a closed circuit. 
At higher initial CO, concentrations, the enrich- 
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COs in initial gas mixture 
Fic. 3. Total reflux operation — 3-component mixture. 


Column 2: AT = 250°C, 2w = 1-80 cm., H,/CO approx. 
constant. 


ment of CO approaches a value one would expect 
for the binary mixture H, CO. The CO, is itself 
concentrated at the bottom of the column. 

The experiments were repeated under total 
reflux conditions with the introduction of CO 
into H, CO, mixtures. (Fig. 4). The maximum 
was again recorded, with an initial CO concentra- 
tion of 22.0% CO in the mixture. The curve 


w 
2) 


4 


€ *le COz in 
(calc. to a basis of CO free gas) 
a 


12) 10 20 30 40 SO 
*fe CO in initial gas 
Fic. 4. Total reflux operation - 3-component mixture. 
Column 2; AT = 250°C, 2w = 1:80 ems. H,/CO, of the 
same order of magnitude (2-56-3-54). 
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did not show such a rapid increase in Ac, i.e. 
improvement in the separation of CO, from H,. 
These results confirm the experimental values 
obtained for the introduction of CO, into H,/CO 
mixtures. 

Consider the result obtained with 15% CO, in 
the initial CO/H,/CO, mixture. Initial con- 
centration of CO on the CO, free gas is 28%. If 
the mixture were a binary of H,/CO, this would 
give an end concentration of 29:1% CO, and an 
enrichment factor of 1-27. In the presence of 
CO,, the CO end concentration, to a basis of 
CO, free gas, becomes 53-0%, with an enrichment 
factor of 2:30. This is a considerable improve- 
ment. 

In the same way, consider the separation of 
CO, from H,. In the absence of CO, the enrich- 
ment would be from 18-5% CO, to an equilibrium 
value of 36-:0% an enrichment fact of 1-58. In 
the presence of CO however, the end concentration 
to a basis of CO free gas becomes 43-8%, with 
an enrichment factor of 2-36. It is again apparent 
that the presence of CO enhances the separation 
of CO, from H, quite considerably. 

It may be concluded from these results that 
any design would have to take into account the 
mutual effect of the molecules, and it cannot be 
accepted that the separation of CO,/H,, CO /CO,, 
and CO/H, mixtures independently would be a 
guide to the behaviour of the three component 
mixture in the column, There is a large density 
change along the column, so that the premises 
used in theoretical work with isotopes would not 
hold for common gases. 

The reasons for the results of these experiments 
on ternary mixtures, are a matter of conjecture 
in the absence of an exact theory for the separa- 
tion of common gases either in binary or ternary 
mixtures. The effect may be due entirely to 
density and viscosity changes improving the 
operating characteristics of the column. This 
opinion would be quite acceptable for the im- 
proved separation of CO from H, by the intro- 
duction of CO,, as the heavier CO, molecules 
improve the thermal syphoning, and the CO,/H, 
system has a better separation factor than 
CO/H,. If however the molecules do not act 


independently, then it would not be surprising 


if the greater ease of separation of CO, from H, 
would affect the CO /H, separation. It does not 
however, and other reasons have to be found. 

A reason which appears to be worthy of some 
consideration, is the possibility that ‘ sweep 
diffusion’ is superimposed on the thermal 
diffusion effect. In this way three processes 
would have to be considered in the design of a 
Clusius type column, (i) thermal diffusion, (ii) 
natural diffusion, (iii) sweep diffusion. Sweep 
diffusion is discussed by CicHeLL1, WEATHERFORD 
and Bowman [17] and is itself a method for the 
separation of gases. In this method a readily 
condensible gas or vapour is introduced into a 
binary mixture. This vapour moves towards a 
cooled surface where it is condensed and removed 
from the system. The movement of the vapour 
has a drag effect on the molecules, the effect 
being different for the components being separated 
The component having the lower natural diffusion 
rate through the vapour is “ swept” out of the 
mixture. 

In the present experiments, it could be con- 
sidered that the CO, being separated more easily 
from H, than the CO, moves rapidly to the cold 
wall, where it is removed from the separating 
zone by the syphon action of the column. In 
this way, superimposed on the movement of the 
CO molecules to the cold wall, is a sweep action. 
The CO will have a lower coefficient of natural 
diffusion through the CO, than through H,, so 
that back-diffusion of CO will be reduced and an 
enhanced separation is obtained over that 
obtainable by thermal diffusion alone. This 
would also explain why the effect is greater for 
the introduction of CO, into H,/CO mixtures, 
than for CO in H,/CO, mixtures. The theoretical 
analysis would be extremely complex, but the 
optimum concentration of the heaviest component 
could be found as above from pilot plant study. 

In the same apparatus results were obtained 
for the effect of withdrawal rate, in continuous 
operation, on the separation of H,/CO/CO, mix- 
tures. The results are given in Table 3, and Fig. 5. 
Under steady conditions of operation, the separa- 
tion is impaired by an increase in withdrawal 
rate. This is shown in Fig. 5, each curve repre- 
senting a different gas mixture. The curves 
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+ (The broken curve represents a binary mixture of H,/CO— 
in the same column under identical conditions) 


|__| | | | 


— 


O 100 200 300 400 500 600 700 800 S00 1000 1100 


Withdrawal rate 7 


Continuous operation — 3-component mixture. Column 2: AT = 250°C, 2w = 1-80 cms, 


intersect the ordinate at the total reflux value 
of the concentration. On the same graph is 
plotted the experimental curve obtained with a 
binary mixture of H, CO with initial concen- 
tration of 15% CO in the same column under 
identical conditions, It is seen that the enhance- 
ment effect of CO, on the separation of H, ‘CO is 
still present as for the total reflux operation. 
Further the values from total reflux operation 
seem to be in agreement with the plotted curves 
from the results of continuous operation. The 
‘maximum ° effect in the separation due to the 
introduction of CO, is still present. The extent 
of the increase in the enrichment over that 
obtained in a straight binary separation, can be 
determined at any withdrawal rate. This increase 
in enrichment is reduced by an increase in the 
withdrawal rate, until at approximately 1,000 
ec/min the separation of H,/CO is as for the 
binary mixture. 

The design of columns for the continuous 
separation of a three component mixture, presents 
an even greater problem than for “ total reflux ” 
as it has been demonstrated in the present 
experimental work, and elsewhere by the authors 
{15]. The characteristic constants H and K are 
influenced by the withdrawal rate for a given 
mixture in a given system under steady state 
conditions. When the other effects associated 
with ternary mixtures are superimposed, the 


problem would probably yield an extremely 
complex theoretical solution, and for design a 
more empirical approach is necessarily based on 
experimental pilot plant data from prototype 
thermal diffusion columns. 


CONCLUSIONS 
(i) The existing and exact and approximate 
theoretical solutions for the separation by 
thermal diffusion of multicomponent mix- 
tures do not apply to common gases, 


The transport equations developed for the 
separation of isotopes, whether binary or 
multicomponent mixtures, cannot be used 
for the design of vertical thermal diffusion 
columns as applied to the separation of 
common gases. 


It would appear that when considering the 
separation of ternary mixtures of common 
gases in vertical columns, the components do 
not act independently. There is a net effect 
governed by the extent of interaction be- 
tween the molecules, and the particular 
physical characteristics of the mixture, such 
as density and viscosity. 


Experiments carried out both under total 
reflux and continuous operation, indicate an 
optimum concentration of a third component 
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for which the separation of the components ence has been made to the mixtures of 
of the mixture is a maximum. Further, the H,/CO/CO,, and the uncertain state of the 
enhancement of the separation of the com- theory necessitates restricting the conclusions 
ponents concerned is a significant one. Refer- to these three gases. 
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The Graetz-Nusselt problem for a power-law non-newtonian fluid* 
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Abstract—The authors show how the Graetz-Nusselt problem in heat transfer theory may be 
extended to non-Newtonian flow. A simple rheological model, * the power law ,”’ is used to 
obtain a partial differential equation for the temperature profiles, to which a semi-analytical 
solution may be found. The temperature profiles are then used to calculate average outlet 
temperatures, as well as Nusselt numbers for several degrees of non-Newtonian behaviour. 
Tables are given for estimating the effect of non-Newtonian flow on the heating and cooling 
of fluids in tubes. 


Résumé— Les auteurs montrent comment étendre le probleme de Graretrz-Nusse.t dans la 
théorie du transfert de chaleur 4 un fluide non newtonien. Un modéle rhéologique simplifié, 
désigné sous le nom de “ loi de puissance * permet d’obtenir, pour la répartition des températures, 
une équation différentielle partielle, qui peut étre résolue semi-analytiquement. Les répartitions 
des températures sont alors utilisées pour calculer une température de sortie moyenne ainsi 
que les nombres de Nussevt pour différents écoulements non-newtoniens. Des tables permettent 
d’estimer l'effet d’écoulements non-newtoniens sur le réchauffage et le refroidissement de fluides 
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circulant dans des tubes. 


l. INTRODUCTION 


Tue solution to the following heat transfer 
problem was presented in 1885 by Graerz [5] and 
again in 1910 by Nussevt [8]: A fluid flows with 
a fully developed laminar parabolic velocity 
profile in the + z-direction in a circular tube of 
radius R. In the region z < 0 the fluid is at a 
uniform temperature 7',. At the section z = 0 the 
fluid passes into a tube of the same radius, the 
walls of which are maintained at constant tem- 
perature T,, greater or smaller than 7',, for all 
z>0. 

Graetz and Nusse.t obtained the temperature 
profiles 7’ (r, z) in the flowing fluid by making the 
following assumptions : (1) steady state has been 
attained, (2) heat conduction in the z-direction is 
negligible in comparison with heat transport in 
the z-direction by the over-all fluid motion, (3) the 
equilibrium physical properties p’ (density) and 
c (heat capacity per unit mass) are independent 


of position, (4) the nonequilibrium properties yu 
(viscosity) and k (thermal conductivity) are 
independent of position, (5) heat produced by 
viscous dissipation is neglected, (6) there are no 
external (body) forces acting on the fluid, and 
(7) the flow obeys Newron’s laws of viscosity 
(which leads to the parabolic velocity profile). 
The calculations of Grarrz and Nusse._t have 
been critically revised and extended by Drew [4] 
and by Yamacarta [9]. 


This paper is concerned with the relaxation of 
the seventh assumption, so that a solution to the 
Graetz-NussE problem for one non-Newtonian 
model is obtained. The non-Newtonian model 
chosen for the calculations is the “ power-law,” 
which has been proven useful for the description 
of pressure drops in flow systems, and for which 
some heat transfer calculations have already been 
made [1]. The flow of many plastics melts may be 
described in terms of this power law. 
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2. ANALYTICAL STATEMENT OF 
THE PropLem 


The steady flow of an incompressible fluid with 
no external forces is described by the set of three 
basic equations of continuity, motion, and 


energy [2] : 


Continuity : 
(V-v)=0 (1) 


Motion : 
p’(v: V)v = — Vp 


Energy : 
pe(v-V)v = —(V-q) —(r: Vv) (3) 


in which v is the local fluid velocity, p is the static 
pressure, T is the temperature, q is the heat flux 
vector, and + is the stress tensor. The term 
(+: Vv) in the energy equation is the heat pro- 
duction due to viscous dissipation. The rdle of 
this term in Newtonian flow in tubes has been 
studied by Brinkman [3] and its importance in 
extrusion of molten plastics has been discussed by 
one of the authors [1]. Under normal conditions 
it is not important and is hence neglected in the 
remainder of this discussion. 

For flow in cylindrical tubes the above equations 


(2) 


in which +,, is the shear stress in the z-direction 
per unit area on an element of fluid surface of 
constant r, and g, is the heat flux in the radial 
direction r. In order to solve these equations we 
substitute for the fluxes r,, and qg, the expressions : 


* Power Law’ : 


Fourter’s Law : 


G = —k — (7) 


or 
in which m and n are constants in the empirical 


Fic. 1. Velocity distribution for the power law non-New- 
tonian fluids with n = 1, 2, 3, @. 


modification of Newron’s law, obtainable from 
analysis of flow data in various types of viscosi- 
meters [7], & is the coefficient of thermal con- 


Table 1. The eigenvalues c,. 


7 


| 


6-582 
39-09 
99-50 

187-9 


7314 | 


1-516 
0-866 
0-634 
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Motion : 
_ dp id 
| 
2 | “61 8635 | 8947 
| 9234 | 74-89 
4 | 215-25 | - | _ 
Table 2. The expansion coefficients B,. 
dv,|" dv, | 
1 | 1-676 1-493 | 3-602 
| | = 
| | | | 
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Table 3. The Eigenfunctions 


1-000 1-000 1-000 1-000 1000) (1-000 1-000 1-000 


0-884 0-766 0-984 0-782, 0-986 0-925 | 0-821 


0/603, 0-646 0-288) 0-938 9-668; 0-269 0-943 0-718 
| 


— — — 


0-234) — 0-315 0-855 0-310 — 0-257) 0-854 0335 — 0218 0-874 0-423 | — 0-095 
i | | | 


0-4 60-738 — 0-110 — 0-392) 0-756 0-045 — 0-403) 0-765 | —0.002 — 0-404 0-782 0-106 — 0-375 


| | 


0-5 0-615) — 0-342) 0-142) 0-638 0-311 — 0-209 0-647 — 0-264 — 0-249 0-670 | — 0-168 — 0-356 


vet. 06 0-483 — 0-432 0170 0-506 |— 0-414 0-118 0517 0-399 0-062) 0543 0-347) — 0-118 


1956/57 O-7 | 0-351 | — 0-398) 0-332 0-374 |— 0-385) 0-310) 0-380 |— 0-400) 0-286) 0-408 0-403) 0-166 


0-8 | 0-224 — 0-285 0-303 0-236 — 0293 0310 0-245 —~ 0301, 0306 0-268 — 0-339 0-299 


09 0107 — 0141 0168 O110) — 0166 O185 — 0153 0-174 — 0-188 0-218 


10 0-000 0-000 0-000 0-000 0-000 000 0-000 o-000 «60-000 0-000 | 0-000 


Table 4(a). Reduced temperature profiles O(p, ¢) for n = 


n=1 n=2 n=83 n = © 
0-0 | 1-000) 1-000) 1-000, 1-000 | 
0-2 0-929 0-380 
03 0-846 
|. 
| | | 
| 
Vavg 
0-10 | 0-20 0-40 0-60 0-80 1-00 
p | 
0-0 0-939 0-701 0-341 0-165 0-079 0-038 0-001 
0-930 0-689 0-335 0-162 0078 0037 0-001 
0-2 0-654 0-317 0-153 0074 0-035 0-001 
0-3 0-845 0-598 0-289 0-139 0-067 0-032 0-001 
0-4 0-765 0-525 0-252 0-121 0058 0-028 0-001 
0-5 | 0-659 0-440 0-210 0-101 0049 0-001 
0-6 0-582 0-347 0-165 0-080 0038 0-001 
0-7 0-395 0-253 | 0-120 0-058 0-028 0-013 0-000 
0-8 0-255 0-162 0-077 0-037 0-018 0-009 0-000 
0-9 0-122 0-077 0-036 0-018 0-008 0-004 0-000 
10 0-000 0-000 | 0-000 0-000 0-000 0-000 0-000 
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Table 4(b). Reduced temperature profiles @(p, ¢) for n = 2. 


Table 4(c). Reduced temperature profiles 0 (p, ¢) for n = 8. 


ductivity. The use of this “power law" has in which p=r/R is the dimensionless radial 
been discussed in detail by Reiner [7]. For variable, and v,,, is the maximum value of v, 
Newtonian liquids m is the coefficient of viscosity at p= 0. Thd velocity profiles for n = 1, 2, 8 
pandn = 1. and o are shown in Fig. 1. Substitution of 

Substitution of eq. (6) into eq. (4) gives an Fourter’s law (eq. 7) and the velocity distribu- 
equation for the velocity distribution which may tion (eq. 8) into eq. (5) then gives the partial 
be integrated to give : differential equation for the temperature profiles : 


= (1 — p)**? (8) in which 6@=(T —T,)/(T, — T,) and 


0-20 0-40 0-60 0-80 1-00 2-00 
0-0 0-926 0-670 0-808 0-140 0-063 0-029 0-001 
0-1 0-923 0-663 | 0-804 0-188 0-063 0-028 0-001 
0-2 0-888 0-629 0-288 0-181 0-059 0-027 0-001 
0-8 0-834 0-577 0-263 0-119 0-054 0-025 0-001 
0-4 0-764 0-513 0-233 0-106 0-048 0-022 0-000 
0-5 0-667 0-485 0-890 0040 0-000 
0-6 0-543 0-346 0-156 0-071 0-082 0-015 0-000 
0-7 0-408 0-256 | O-115 0-052 0-024 0-01! 0-000 
0-8 0-262 0-162 0-073 0-083 0-015 0-007 0-000 
0-9 0-124 0-076 0-084 0-015 0-007 0-003 | 0-000 
1-0 0-000 0-000 0-000 0-000 0-000 0000 (0-000 
Vol 
6 
195¢ 
en | con 0-60 0-80 1-00 2-00 
| | 
0-923 0-651 0-285 0-124 0-054 0-023 0-000 
0-014 0-641 0-281 0-122 0-053 0-023 0-000 
0-886 0-612 0-268 0-116 0-050 0-022 0-000 
0-827 0-559 0-244 0-106 0-046 0-020 0-000 
0-764 0-508 0-218 0-095 0-041 0-018 0-000 
0-666 0-427 0-185 0-080 0-035 0-015 0-000 
0-547 0-848 0-148 0-064 0-028 0-012 0-000 
0-411 0-258 0-108 0-047 0-020 0-009 0-000 
0-268 0-163 0-070 0-080 0-013 0-006 0-000 
0-129 0-078 0-088 0-015 0-006 0-003 0-000 
0-000 0-000 0-000 0-000 0-000 0-000 0-000 
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Table 4(d). Reduced temperature profiles 0 (p, 0)for n = oo. 


ae 


| 


| | 


Table 5. Average outlet temperature 


| | | | 
1 2 3 x 


"max 


0-579 
0-395 
0-190 
0-091 
0-044 
0-021 
0-001 


0-558 
0-370 
0-168 
0-076 
0-035 
0-016 
0-000 


Table 6. Nusselt numbers. 


| 
| 


l = (k/p’cv,,,,R*)z. Eq. (9) is to be solved for 
the boundary conditions: @=1 at (=0; 
@=0forp=1; and = 0 atp = 90. 


8. SOLUTION FOR THE 
TEMPERATURE PROFILES 


Equation (9) is solvable by the method of separa- 
tion of variables, and the radial portion of the 
problem is of the Sturm-Liouville type. Hence 
the temperature profiles must be of the form : 


in which the factor (— 1)'*! is included so that 
the tabulated B; are positive. The ¢,(p) are the 
eigenfunctions (orthogonal with respect to the 
weight function p(1 — p"*') in the range p = 0 
to p = 1) obtained from the solution of : 


dp 
the c; being the eigenvalues. The ¥,(¢) are 


= exp — ef 


and the B; are determined from the condition at 
= 0 to be: 


(12) 


1 
faa —p"*") pdp 
B,=° 


—p"*")pdp 
0 


time 
\ Pave 
\ 0-10 0-20 0-40 0-60 . 1-00 2-000 
p | 
| 
0-0 0848 0502 159 0-050 0016 0-005 0-000 
O1 0-839 | 0-494 0-156 0-049 0016 0-005 0-000 
0-2 0-811 0473 149 0-047 0015 0-005 | 0-000 
0-3 0764 0-439 | 0-189 0-044 0-014 | 0-004 0-000 
0607 0-039 0012 0-004 0-000 
05 0-610 0-837 (0-106 0-033 0011 0-000 
0-6 0506 0275 0-086 0-027 0-009 0-000 
0-7 0-206 0-065 0-020 0-006 0-002 0-000 
0-8 0258 0-136 | 0-042 0-013 | 0-004 | 0-001 0-000 
0-9 0127 0066 0-021 0-006 0-002 0-001 0-000 
10 0-000 | 0-000 (0-000 0-000 0-000 
| 
| | | 
ol 
o2 | 0-351 O-216 
06 | | 0-066 0-021 
| 0-020 0-007 
10 | 0-012 | 0-002 
n 1 2 3 | re 
o1 | | «02 | #50 | 606 
| 303 | 315 | 325 | 892 = 
04 203 908 2-33 
0-6 152 | | 1-56 1-63 
10 0-98 0-98 0-99 1-00 
2-0 0-50 | 0:50 0-50 0:50 
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Hence the problem is reduced to the deter- 
mination of the eigenvalues c; and the eigen- 
functions ¢,(p) for several values of n of interest. 
For x = 1 the most extensive calculations are 
those of Yamacata [9], for n = 2 calculations 
made earlier by one of the authors may be used [1], 
for n = 3 new calculations [6] are reported here ; 
for n = o (flat velocity profile) the ¢; are just 
the zero-order Bessel functions J, (\/cip). The 
¢; may be obtained as series solutions : 


= 


in which the 6; 4 are arbitrarily chosen to be 
unity and the ‘higher coefficients are obtained 
from the recursion formula : 

= (0, — 54, 5-2) (15) 
The values of the ¢; are obtained from eqs. (14) 
and (15) by use of the boundary condition at 


(14) 


=1 which requires that = 0. 
0 


The values of the ¢; are een in Table 1, and 
the corresponding values of B; may be found in 
Table 2. The eigenfunctions and reduced tem- 
peratures are presented in Tables 3 and 4 respec- 
tively. The tabulated values of ¢, and B; for 
i = 8 and n = 2,3 are probably not so accurate 
as the number of figures would indicate, inasmuch 
as the accuracy was limited by the use of a 10-bank 
desk calculator. The tabulated values of the 
reduced temperatures are good to about 2°,. Some 
sample temperature profiles are shown in Fig. 2. 
From the tables and graphs of reduced tem- 
perature profiles one may estimate the tempera- 
ture distribution for the non-Newtonian fluid 
obeying the power-law model. 


4. AveRAGE OuTLET TEMPERATURE 
AND NusseELT NUMBERS 

For many problems the average outlet tempera- 
ture of the fluid after passing through a length 
of pipe from z = 0 to z = z may be desired. This 
temperature is defined by : 

271 

Jo £) o(p) pdpd¢ 


= 


(16) 


"J 


10 


Fic. 2. Reduced temperature @(p, () for n = 1, 2, 3, @ 
for (Umax avg) = 0-6. 


Values of 6,,, are given in Table 5 as a function 
of the reduced length { multiplied by (0,.4./ayg)s 
which ratio has the values 2, 5 3, 3/2, 1 forn = 1, 
2, 3, «© respectively. This table allows a com- 
parison of the reduced temperature change for 
fluid with the same thermal diffusivity and the 
same average velocity but with different degrees 
of deviation from NewrTon’s law of viscosity. 

The results may also be summarized in terms 
of the Nusselt number for heat transfer. First a 
heat transfer coetticient h is defined by the 
relation : 


q =h-2eRz-(T, — T,) 


in which q is the heat flow per unit time through 
the cylindrical surface of area 27Rz. The heat 
flow through this surface must by an over-all 
energy balance be the same as the difference 
between the sensible heat input by flow at z = 0 
and the output at z =z. That is: 


q = (Ty — (2)) (18) 


Combination of (17) and (18) along with the 
definition of the Nusselt number as Nu = 2Rh/k 


gives : 
1— 
Uniex/Vavg 


The Nusselt numbers so calculated are shown in 
Table 6 for n = 1, 2,8, 0. The values of Cun 


(17) 


(19) 
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and Nu were obtained by numerical integration NOMENCLATURE 
with Simpson’s rule. bj, ; = coeflicients in (14) 
c = heat capacity per unit mass 

; = eigenvalues in (11) 
= coeflicient of thermal conductivity 
= constants in power law model 
= static pressure 
= heat flow in (17) 
= heat flux, radial component (5, 7) 
= heat flux vector 
= radial co-ordinate 

velocity vector 

z-component of velocity 


5. ConcLUDING REMARKS 


The foregoing calculations represent only the 
first step in the study of heat conduction in non- 
Newtonian flow systems. Similar calculations 
need to be made for other non-Newtonian flow 
models (such as the Bingham plastic model, for 
example). Furthermore, the temperature depend- 
ence of the viscosity needs to be taken into 


<1 


account. Also, the tube-flow problem here dis- average flow velocity 
cussed is but one of many geometrical arrange- max = ™aximum flow velocity 
ments which deserve attention. = axial co-ordinate 
At the present time it is not possible to make 
any comparisons with experimental data. There = padius of tube 
VOL. seem to be no data on heat transfer in laminar = temperature 


non-Newtonian flow, although such systems are = wall temperature 
of considerable industrial importance. There is = inlet temperature 
also urgent need for experimental data on thermal yg ™ average exist temperature 


= dimensionless axial co-ordinate 
conductivities and heat capacities of non-New- = dimensionless temperature 


tonian fluids. = dimensionless average exit temperature 


= N . 
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Abstract — Measurements have been performed on the rate of solution of rather short surfaces 
of benzoic acid in a water film in laminar and turbulent flow, For the experimental conditions 
this process may be treated as the diffusion of solute from a plane surface into a laminar liquid 
flow with a constant velocity gradient. A reasonable agreement is found between the velocity 
gradients derived from the measurements and those which follow from otherwise known charac- 
tersitics of falling films. 


Résumé—-On a mesuré la vitesse de dissolution dans un film d'eau, de plaquettes d'acide 
benzoique. On a opéré dans les domaines d'écoulement laminaire et turbulent. 
Dans les conditions expérimentales réalisées, ce probléme peut étre traité comme la diffusion 
d'un solute a partir d'un plan a travers une couche laminaire a gradient de vitesse constant. 
Les valeurs du gradient de vitesse calculées a partir de ces expériences sont en bon accord 


avec celles que l'on peut déduire des autres caracteristiques des films tombants. 


1. INTRODUCTION 


Tuts study was undertaken because of interest in 
the resistance to heat or mass transfer across a 
liquid film flowing over a flat surface. In parti- 
cular in the turbulent region such information 
might be of use for predicting the rate of heat 
transfer between a vapour and a cooled surface 
at high rates of the condensate flow. In principle, 
an experimental determination of the total 
resistance of the liquid film would require the 
establishment of a constant driving force over its 
two boundaries and the measurement of the 
resulting heat or mass flux at a point where the 
final temperature or concentration distribution 
has been reached. 

In the present investigation a different approach 
to the experimental problem has been taken, and 
as a consequence, the information originally aimed 
at was not obtained. The experimental procedure 
was based on the assumption that the main 
resistance for transfer through a turbulent film 
would be located near the solid surface. This is 
not unreasonable since Brérz [1] has shown with 
experiments on the spreading of dyes and on gas 
absorption in turbulent films that there exists a 
fairly high coefficient of mixing or “ eddy diffu- 


sivity near the free surface. Thus, the experi- 
ments reported here had the purpose of deter- 
mining the mass transfer rate between the solid 
surface and the bulk of the liquid film. 


To this end a falling film of water was made on 
an inclined flat metal plate. At a downstream 
position, where the hydrodynamic flow was 
sufficiently established, a part of the plate was 
substituted by a soluble solid (benzoic acid) and 
its rate of dissolution was measured. The length 
of the soluble surface was varied between 5 and 
80mm. Now, owing to the small diffusivity of 
benzoic acid in water the depth of penetration of 
the solute into the liquid under these circumstances 
is so small that only the conditions in the laminar 
boundary layer are of influence on the rate of 
mass transfer. Thus these experiments only give 
information on the nonsteady diffusion into a 
laminar liquid flow with a constant velocity 
gradient and the further discussion will be res- 
tricted to this particular case. For obtaining data 
on the behaviour of the liquid beyond this 
boundary layer, much greater lengths of soluble 
substance would have been required. Experi- 
ments of this kind have recently been reported 
by and Hurt [2]. 
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2. TnEory 


The problem of heat penetration from a surface 
into the laminar boundary layer of established 
flow past the surface has already been solved by 
Exsex [3]. His result can equally well be used for 
the analogous case of diffusion. Since, however, 
Exser’s derivation is rather elaborate and a more 
elegant way of solving this problem exists, the 
latter will be given below. 

The following problem, then, has to be solved. 
A flat plate of arbitrary breadth consists for 
x < 0 of nonsoluble material and for + > 0 of a 
substance with solubility c* (Fig. 1). A liquid 
with solute concentration ¢ = 0 is flowing past 
the surface in the positive 2 direction. The flow 
is laminar and it has a constant velocity gradient 
dv/dy = a and a velocity v = 0 at y= 0. The 
average rate of mass transfer is to be found be- 
tween 0 and 2. 


Fic. 1. Conditions for which eq. (1) is solved. 


Neglecting the molecular diffusion in the 
direction of flow, we have the differential equa- 
tion : 

ve de 


D (1) 


with the boundary conditions : 


c= 0 forr = Oandy > 0, (2) 
c= c*forz >Oandy = 0. (3) 
Equation (1) has a special solution of the source 
type : 
= exp (— y'/2), (4) 
where 
9D 
a 


S is a constant which represents the source 


strength. 
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Solution (4) already obeys the boundary con- 
dition (2). By a superposition of infinitely small 
sources dS({) a general solution is obtained of the 


form : 
= 


¢= (z — dS(Z). (5) 


The function dS({) follows from boundary con- 


dition (3) : 
2/3 ds d 6 
on (2 — a, (6) 


The latter equation yields : 


dS 


(1) 


Substitution of (7) into (5) gives for the con- 
centration distribution : 


~ 
— 


c* 
| 
3 3 
Since we are interested in the rate of dissolution 
of the solid, the concentration gradient at y = 0 
has to be calculated from (8). This gives : 
3 
3 


The value of the mass transfer coellicient, k, 
averaged over the distance 2 then becomes : 


ct (9) 


z 


0-808 (=) (10) 


Exactly the same result can be derived from 
Exser’s solution of the problem. 

A further development of the very general 
equation (10) depends on the hydrodynamic 
circumstances which indicate in what terms the 
velocity gradient a has to be expressed. For 


0 
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6 
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x 
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gravity flow of a liquid film along a flat surface 
we have in the steady state : 


(**) (11) 
dy/y-o0 


where g, is the component of the acceleration by 
gravity in the direction of flow x. The film 
thickness 5 may be expressed in terms of the 
liquid flow rate I). For laminar film flow we 
have the relationship : 


This relation may also be written as : 


According to many authors (12) or (12a) may be 
used up to Re — about 2,000, although already 
at low Re numbers rippling occurs and the velocity 
profile deviates from the theoretical one. 

From experiments on turbulent films falling 
along a vertical wall Brirz [1] has shown that the 
film thickness may be described by the relation- 
ship : 


2\ 1/3 
v 
in which the viscosity no more occurs as a 
variable. Eq. (13) may also be written : 
1/3 2/3 
v | 
Some controversy seems to exist in this field. 
Dux er and [4] report that already at 
Re = 1,000, 8 becomes greater than the value 
predicted by the laminar theory. However, 
Jackson [5] recently showed that for a vertical 
water film eq. (12) applies at least up to Re — 4000, 
In view of these uncertainties the equation for the 


film thickness in laminar and turbulent flow along 
a flat plate will be written here : 


* A similar result has been reported by Kamer and 


(14) 
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where F is assumed to depend on Re only and 
not on the slope of the plate. 

Substitution of (11) and (14) into equation (10) 
gives the final expression for the mass transfer 
coefficient k, averaged over a distance 2, between 
a flat surface and established laminar or turbulent 
film flow : 


For the experimental verification of this 
expression, it is essential that the penetration of 
solute does not extend beyond the region where 
the velocity gradient is constant. 


(15) 


1/9 
) = 0-808 F', 


4. EXPERIMENTAL 


The experiments were performed in a flat open duct which 
was 100 mm wide (Fig. 2). Its slope could be adjusted and 
for these measurements three values were chosen of the 
angle, a, between the surface and a horizontal plane : 80, 
45 and 6 degrees. At a distance of 330 mm downstream 
slabs of soluble material could be inserted in the manner 
indicated in Fig. 2. For these experiments benzoic acid 
was used and the slabs were cast separately, 100 mm wide 
and with variable length. Great care was taken that the 
slabs had a flat surface and that this surface could be fitted 
flush with the bottom of the duct. Before use the slabs 
were painted with a very thin layer of plastic paint (quick 
drying and smooth) except for a rectangular surface on the 
top side which was 80 mm wide and of which the length 
2 was varied between 5 and 80 mm. The rate of dissolution 
of the slabs was measured at four different water flow rates 
which were measured with a calibrated rotameter. 


Woter iniet 


. 


Paint 
/ 


L 90mm" 


Fic. 2. Apparatus for the measurement of mass transfer 
to film flow, 


(12a) 
6 
z 
surface 
— 
|| 
Spring 
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The average mass transfer coeflicient k over the length 
av, was calculated with the formula 


AM 
ba: (c* — Cen) : 


(16) 


the area of the soluble surface and ¢ the 
time of exposure. The difference in mass A.M of the slabs 


where br 


before and after the experiment was obtained by weighing. 
Before weighing the slabs were dried during half an hour at 
orc, For the dissolution time a compromise had to be 
chosen. Too short a time would result in inaccurate AM 
values and at too long times too much of the surface would 
be dissolved. According to circumstances ¢ varied between 
3 and 15 minutes. Since for the x values used the degree 
of saturation of the water was extremely small, ¢,,, could 
be neglected with respect to the solubility of benzoic acid 
in water, c*. The latter value, as well as the diffusivity of 
benzoic acid in water at the experiment’s temperature, was 
taken from the paper by Linron and Surrwoop 

The experimental results have been collected in the 
Tables 1, 2, and 3 for « 
In the right-hand column the dimensionless group con- 
taining & is shown. According to eq. (15) this group should 
depend on ReyNoups’s number only, so it should be con- 
stant for a series of experiments where « and I), are kept 
constant. For most combinations chosen this appears to 
be the case within a standard deviation o of 5 to 10°. 
Only for « = 80 deg. and the two highest flow rates 
(Table 1) this group shows a significant trend as 2 is varied. 
From 14 duplicate runs (marked * in Tables 1 and 2) it 
was found that & had been measured with a standard 
deviation of about 8°. 

As a summary of the results Table 4 shows the quantity 
F1/3 which was calculated from the mean values in Tables 
1 to 3 according to eq. (15). Considering the standard 
deviations, there is no significant dependence of F!/8 on 
a. There is, on the other hand, a significant correlation 
between and Reynoups’s number, 


80, 45 and 6 degrees respectively. 


5. Discussion 
Since the accuracy of the experimental k-values 
is not great and since k does not strongly depend 
on the velocity gradient (compare eq. 10) it is 
hardly possible to use these measurements for 
predicting the thickness of falling films. On the 
other hand, it is of interest to compare the values 
found for F'* with those which can be deduced 
from the various existing correlations. Such a 
comparison is presented in Fig. 3. On the whole 
there is a fairly good agreement between the 
experimental values and those from the laminar 
theory (Re < 2,000), and from the correlation (13) 
for turbulent flow. Still, one might want some 
explanation as to why the experimental values of 


45 


I'* at Re ~ 3,500 and 5.000 are low and those 
at Re ~ 7,000 high. 


4 


a =80 
a=45 
a=6 
Laminar theory, eq{12) 
Dukler [4] 
Brotz [1] 

| | 
ie) 2000 4000 6000 8000 


Fic. 3. Comparison between values of F'3 (Table 4) and 
those derived from measurements of the film thickness. 


Table 1. Mass transfer coefficients k, for « = 80°; 
T = 21°C. 


v 


k vas 


3-13 

2-97 

294 

2-60 

2-53 | 

2-12 | mean = 

1-73 

2-62 

2-22 

2-15 

200 

2-02 | 

1-74 mean = 2-06 
1-67 = 0:30 


194 
1-90 
1-91 
1-84 | 

167 mean = 1-90 
2-14 a=0-19 


| 


Note : k value marked * is an average of two measurements 


| 
| 2 
3 
| 
4 
71 80 478 | 
| 0 500 
0 | 565° 
30 | 5-50* 
20 613% 
| | | 
5 | 665 | 
60 357° | 
30 $23 
20 490 | 
10 CO 
| | | 
206 
| 6 | 320 
| | 866 
| | | 
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Table 2. Mass transfer coefficients k, for x = 45°. 


zx 
10°? m 


| 
| 


ows 


| 


k 
10°* m?/sec 10 5 m/sec | 
17-1 204 
3-47 2-22 
+76 2-56 
3-54 2-06 | 
| 468 2-34 
448 2-27 
6-12 2-63 | 
5-73 | 2-30 | 
6-45 2-20 | 
7-01 2-25 | mean = 2:37 
| 8-62 2-35 o= 024 
12:8 2-27 1-83 
| 3-18 2-20 
3-34 2-08 
3-07 1-96 
| 20 170 vol 
3-40 1-98 
4-00 1-99 195€ 
| 3-80 1-93 
| 5-25 224 
+70 1-89 
5-15 1-75 
5-87 1-88 | mean = 1-07 
8-02 2-16 o = 0-16 
8-3 2-22 1-79 | 
2-45 1:53 
2-78 | 1-58 
2-79° 1-79 
2-97 1-59 
3-25 1-89 | 
3-54 1-60 
3-30° | 1-67 
424 181 | 
+90 1-96 
463 | 1-56 
5-24 | 1-68 | mean = 1-70 
6-10 1-64 | = 0-17 
ae 2-24 1-85 
2-56° 1-79 | 
2-81 1-77 
2:46 1-54 
38-07 1-65 
2-74° 1-59 | 
3-27 | 1-63 
8-74° 1-91 
3-54 1-55 
441° 1-78 
| 4-83 | 1-64 
5-28 1-69 mean = 1-70 
6-37 1-71 = 0-12 
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Table 3. Mass transfer coefficients k, for « = 6°; 
T 
104m? /sec m 10-5 m/sec | £,° D® 
12-8 80 2-16 2-34 
wo 207 
2-32 2-00 
20 286 1-96 
10 | 3-47 187 mean 2404 
5 4-63 200 o = O15 
8-3 1-89 204 
60 211 207 
222 1-00 
20 2-78 1-90 
10 3-33 180 
5 3-04 1-69 o = 0-13 
80 165 «1-79 
| 6o 198 14 
| 0 2-05 1-76 
| | 247 1-69 
| 10 3-28 | 178 mean = 1-79 
re. 416 1-79 o = 0-08 


Table 4. calculated from experimental results. 


2 F! 3 
v deg. 

7000 80 3-2 4 06 
6500 45 340 4 0-4 

5200 80 26 404 

4900 45 24402 

5200 6 25402 

3400 80 23402 

3200 45 21402 

3400 6 2-25 4 O15 
1500 45 21401 

1600 6 22401 


A reason for the observed discrepancy might be 
that under the experimental circumstances the 
entrance length of 0-33 m was too small for a 
sufficient establishment of the final flow pattern. 
This would result in too small F"* values. It ean 
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be shown that the characteristic length for estab- 
lishing the final flow distribution is proportional 
to 1?" g'*, Thus the effect, if present, would 
be most pronounced at the higher Re numbers 
and a small value of x. Since the results do not 
show such a discrimination, it is thought that this 
effect was of no great importance. 

Too high values of k, and thus of F'*, might 
be obtained if the solute penetrates deeper into 
the liquid than the distance over which a laminar 
boundary layer with a constant velocity gradient 
exists. From eq. (5) it may be deduced that the 
depth of penetration is proportional to (Da/a)'/*. 
On the assumption of von KArmAn’s velocity 
profile the distance over which the boundary layer 
is still laminar would be about 5 (va)*. Thus at 
high values of 2 and a such an effect would be 
found first. This might partly explain the less 
satisfactory results presented in Table 1. Although 
in the measurements at high values of 2 reported 
there the solute may not have reached the limit 
of the “laminar” boundary layer, it is well 
possible that it has penetrated into that layer to 
such an extent that the diffusion is already 
increased by small eddies whereas the velocity 
gradient is still constant. 
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NOTATION 


a = dv/dy (see™!) 
¢ = concentration (kg m~*) 
D = molecular diffusivity (m? sec~!) 


F = function of Re, eq. (14) 
g, = acceleration by gravity along «(= g sin (m sec™*) 


k = mass transfer coefficient (m see~!) 

T = temperature (Cc) 
v = velocity (m sec™!) 
2 = coordinate in flow direction (m) 
y = distance from wall (m) 
a = angle between z-direction and horizontal (deg.) 
I’ = gamma function 

I, = volumetric flow rate (m3 sec™! m=!) 
8 = film thickness (m) 
v = kinematic viscosity (m* see!) 

Re = 


(1) 
(2) 
(3) 
(4) 
[5) 
(6) 
(7) 
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